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Summary 

A high-fidelity blade-element mathematical model for the AH-64A Apache Advanced Attack 
Helicopter has been developed by the Aeroflightdynamics Directorate of the U.S. Army’s Aviation and 
Troop Command (ATCOM) at Ames Research Center. The model is based on the McDonnell Douglas 
Helicopter Systems’ (MDHS) Fly Real Time (FLYRT) model of the AH-64A (acquired under contract) 
which was modified in-house and augmented with a blade-element-type main-rotor module. This report 
describes, in detail, the development of the rotor module, and presents some results of an extensive 
validation effort. 


Introduction 

High-fidelity simulation models of helicopters are needed for a variety of tasks including han- 
dling qualities evaluations, pilot training, simulation of life-cycle upgrades, and accident investigations. 
Many different approaches are currently used to develop mathematical models for helicopters. These 
approaches can be organized into three distinct categories (ref. 1): 1) analytical models, 2) identified 
models, and 3) combinations of 1 and 2. Analytical models rely on dynamic and aerodynamic theories 
and usually attempt to model each component of the helicopter individually. The combined, end-to- 
end, responses of such models are generally accurate for the dominant responses and nonlinearities, but 
are often imprecise in modeling rotor dynamics, especially in the off-axis (ref. 2). Identified models, 
on the other hand, use available vehicle- response data (collected through flight testing) to generate 
models which accurately characterize the end-to-end responses of an existing aircraft. Such models are 
commonly used during prototype testing for optimization of flight control systems, where very accurate 
response prediction is critical. In addition, they are very useful bases for simulation model validation, 
as is demonstrated herein. Identified models generally do not attempt to treat the components of 
the helicopter individually and their region of validity is restricted to the configuration and the linear 
response at the flight condition of the identification. Further, there is no way to identify models of 
aircraft still under development (prior to first flight) since there can be no access to flight-test data for 
such vehicles. 

Component type models are, therefore, the only way to simulate helicopters during the design 
phase, and provide the best way to simulate aircraft over their entire flight envelope. The latter is 



especially true at extreme conditions, such as may be encountered in accident investigations. Aside 
from their applicability to the entire flight envelope, however, there is a further advantage to making the 
investment of time and resources needed to develop a component type model of an existing helicopter. 
This advantage is the availability of flight-test data, which allows engineers to perform the extensive 
testing and comparisons of simulation and aircraft responses that are necessary to develop new modeling 
methods or enhance existing ones. These new or enhanced methods can then be applied to the modeling 
of both existing and future helicopters, improving our modeling capability overall. 

In the late 1980s the U.S. Army initiated efforts to obtain a high-fidelity component type model 
for the AH-64A Apache helicopter (fig. 1). The effort was prompted by the success of two pilot-in- 
the-loop UH-60 accident-investigation-simulations conducted at NASA Ames. Those simulations were 
conducted using a component-type model of the Black Hawk known as GenHel, developed by Sikorsky 
and expanded/improved at Ames (ref. 3). A model with comparable complexity and fidelity was sought 
for the AH- 64 A. 

As a result of a competitive procurement effort, the best available model for the AH-64A, the Mc- 
Donnell Douglas Helicopter Systems’ (MDHS) Fly Real Time (FLYRT) (ref. 4), was obtained including 
validation-flight-data in hover and forward flight. FLYRT reponses were shown to match the provided 
flight data reasonably well. However, a great shortcoming of FLYRT was perceived to be its main-rotor 
module which is based on a map- type approach. This modeling approach was originally conceived in 
the late 70s to allow real-time operation with the limited computational capabilities available at that 
time (ref. 5). The approach relies on a pregenerated map to determine the rotor’s quasi-steady response 
at each cycle. First-order-lag approximations to rotor dynamics are then added to these quasi-steady 
values to complete the rotor response. The problem with this technique is one of complete dependency 
on and restriction by this pregenerated map, especially since the code for the generation of more maps 
was not provided as part of FLYRT. The restriction was deemed surmountable since the computational 
power of current computers makes it possible to run blade-element models in real-time, without the 
need for the generation of a map. 

A blade-element-type main-rotor module was therefore developed to replace the map-type module 
in FLYRT. The goal was to achieve the flexibility that a blade-element-type module, unburdened by 
the restrictions of a pregenerated map, would provide. Further, the access that a blade-element-type 
rotor module allows to the actual physical parameters of the rotor would make it possible to introduce 
corrections, enhancements, and new theories as they are developed. It would therefore allow for con- 
tinual improvement, which would not have been possible with the map- type rotor module. The new 
rotor was incorporated into FLYRT, along with additional modifications as necessary. The new model 
is known as the Blade-Element Model for APache (BEMAP). BEMAP has been validated in compar- 
ison with available flight data in hover and forward flight in both time and frequency domains. The 
frequency-domain validation methods draw heavily on the tools contained in CIFER fi (Comprehensive 
Identification from FrEquency Responses) (ref. 6), as will be discussed. This report describes, in detail, 
the derivation of the rotor equations for the blade-element rotor module. It also presents some results 
from the validation effort. In every case, FLYRT responses have been included for comparison. Note, 
however, that the intent is not to show superiority of one model to the other but to show that the 
blade-element rotor module has been derived and implemented correctly. 
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Mathematical Model Description 

To facilitate the introduction of the blade-element rotor module, FLYRT was first extensively modi- 
fied. The modifications included: 1) reorganization of the code to improve modularity, 2) improvements 
in its input/output capability, and 3) addition of a plotting capability not provided under the contract. 
The basic structure of the new rotor module is based on the main-rotor module of GenHel (ref. 7), 
with the relevant equations for modeling the AH-64A derived with the aid of the symbolic manipula- 
tion program MACSYMA (ref. 8). The derivation of the rotor equations was part of an earlier effort, 
initiated by Chen (ref. 9), to develop complete blade dynamic equations for fully articulated configu- 
rations which use one of the three hinge arrangements most commonly used, i.e., lag- flap-pitch (1-f-p), 
flap-lag-pitch (f-l-p), and flap-pitch-lag (f-p-1). Chen derived the inertial portion of the equations using 
a Lagrangian approach, making no high order or small angle assumptions. MACSYMA was later used 
to rederive the same equations using a Newtonian approach, again making no high order or small angle 
assumptions. The two sets of equations were compared and were found to be in general agreement. 
Though the Apache uses an f-p-1 hinge sequence, the simpler f-l-p sequence was used in the new rotor 
module. The choice was made in order to avoid the added complexity that would have resulted from 
treating blade pitch as a degree of freedom. The reasoning was that given the small lead-lag angles, 
the increased accuracy achieved by using the actual f-p-1 sequence would be minimal. 

The new rotor was then integrated into FLYRT to create BEMAP. This also required the re- 
placement of the trim module and the modification of the equations-of-motion (EOFM) module, as 
will be discussed later. The modules representing other components of the Apache helicopter, i.e., 
fuselage/empennage/ wings, vertical tail/tail rotor, horizontal stabilator, and landing gears were used 
directly from FLYRT with minimal changes required for implementation. The Digital Automatic Stabi- 
lization Equipment (DASE) model used in FLYRT was also retained, however, extensive modifications 
to allow introduction of control inputs from flight-data were made. BEMAP is, therefore, a version of 
FLYRT which has been: 1) equipped with a blade-element main-rotor module, 2) upgraded with new 
trim and modified equations-of-motion modules, 3) restructured to improve modularity, 4) enhanced 
with flight-data-control-input access capability for flight-data comparisons, and 5) updated with a 
versatile plotting option. 

The development of the new blade-element model will be discussed in detail in the next section. 
The trim module and the modifications to the equations-of-motion module will also be briefly discussed 
subsequently. No attempt is made in this report, however, to redocument the modules which remain 
functionally unchanged or only slightly changed from FLYRT. The description and documentation of 
those modules can be found in reference 4, generated by MDBS as part of the contract to deliver 
FLYRT to the Army. 


Development of the Rotor Module 

In order to minimize the development time of the blade-element main-rotor module, the main- 
rotor module of Sikorsky/ Ames GenHel was used as a basic structure. To develop an AH-64A version 
of this rotor, it was necessary to: 1) modify the module to allow noncollocated flap and lead-lag 
hinges, 2) replace all UH-60-specific data with AH-64A-specific data, 3) incorporate new aerodynamic 
coefficient tables and new table-lookup and interpolation routines, and 4) replace all UH-60-specific 
equations with AH-64A-specific equivalents. 
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Like the GenHel main- rotor (ref. 3), the BEMAP rotor considers flapping, lead-lag, rotor rotational 
speed, and inflow degrees of freedom. An equal-annuli method (ref. 7), modified for noncollocated 
hinges, was implemented to divide each blade into elements. Equations for the local velocity and 
acceleration of each element were then derived based on aircraft and blade (lead-lag, flapping, and 
rotational) motions and the element’s position with respect to the aircraft Center of Gravity (C.G.). 

The local velocity, along with local inflow and wind, determine the local angle of attack and Mach 
number for entry into the aerodynamic coefficient tables. The coefficients are then used to calculate 
components of the aerodynamic force and moment per blade element. Summing these over all the 
blade- elements results in the aerodynamic forces and moments per blade. No dynamic twisting or 
bending of the blades were modeled. The preformed linear twist of the blades, however, is represented 
through adjustments of the blade-element pitch angle. 

The local acceleration at the blade element was used to derive equations for the inertial forces and 
moments per blade. Unlike the case of aerodynamic forces, the blade elements used for the inertial 
derivations were differential elements and analytical integration, rather than numerical summation, was 
used along the blade span. No simplifying assumptions were made in the derivation of the equations 
(aside from the assumption of rigid blades) and, therefore, the equations are quite complex. They, 
however, provide a good tool for exploring the effects of higher order terms (usually dropped) on the 
fidelity of the rotor module. 

The derivation of the equations for the rotor forces, rotor moments, and the coupled flapping 
and lead-lag equations of motion are provided in the following subsections and appendix A. Simple 
MACSYMA macros developed to aid in the derivations have also been provided (app. B). The aerody- 
namic portion of the equations follow the methods used in the GenHel main rotor (ref. 7), modifying 
for the AH-64A as necessary. Following the GenHel structure, the BEMAP rotor module contains 
its own integration algorithm which gives flapping and lagging positions and velocities in the rotating 
frame. The AH-64A specific data and the aerodynamic-coefficient tables were obtained from the “Air 
Vehicle Technical Description Data for the AH-64A Advanced Attack Helicopter,” reference 10. For 
angles of attack between -5 deg and +30 deg the values of section lift and drag coefficients are from 
wind-tunnel tests and are functions of local angle of attack and Mach number. For angles of attack from 
+30 deg to +355 deg the values of section lift and drag coefficients are from C-81 equations (ref. 11), 
are independent of Mach number, and the lift coefficient data include the increase in maximum lift due 
to dynamic stall. 

The inflow components are calculated using the Pitt/Peters inflow model implemented as described 
by Peters and HaQuang in their AHS Technical Note (ref. 12). The model, which is based on unsteady 
actuator-disk theory, is valid for forward flight as well as hover and uses coefficients of aerodynamic 
thrust, pitching moment, and rolling moment to calculate the three induced velocity states. 

Derivation of Rotor Equations 

As mentioned previously, the symbolic manipulation program MACSYMA was used to derive the 
rotor equations using a Newtonian approach. Figure 2(a) shows the coordinate systems used in the 
derivations. As can be seen, the model assumes a fully articulated rotor hub using flapping and lead-lag 
hinges in a flap- lag-pitch hinge arrangement, the flapping hinge being closest to the center of rotation. 
From the rotor shaft up, the rotating-shaft (rs) frame has its origin at the rotor hub and its y-axis 
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along the blade segment between the hub and the flapping hinge and rotating with the rotor. Frame 3 
is similar to the rotating-shaft system except that its origin is at the flapping hinge. Frame 2 also has 
its origin at the flapping hinge but has its y-axis along the blade-segment between the flapping and 
the lead-lag hinges, i.e., rotated about the negative x-axis through (3 (flapping up is positive). Finally, 
frame 1 is similar to frame 2 except that its origin is located on the lead-lag hinge. 

Velocity and acceleration vectors at the rotor hub— The translational and rotational velocity 
and acceleration vectors at the main-rotor hub due to the motion of the helicopter are independent of 
the hinge sequence used for the rotor. These vectors are the result of the translational and rotational 
velocity and acceleration of the aircraft’s C.G. and its separation from the rotor hub (refer to fig. 2(b)). 


Let: 


( \ 

u 


Vr 


cgb 


W 

v / 


> = Inertial translational velocity of aircraft C.G. (body axes) 


U), 


cgb 


V 

q } = Inertial rotational velocity of aircraft C.G. (body axes) 
r 


Xh 

l h b = { Vh ^ — Hub location relative to C.G. (body axes) 

The translational velocity at the rotor hub, in body axes, may then be calculated from: 


Vhb ~ ^cgb T bJ C gb ^ Gfc 


substituting, we get: 


(i) 


v hb ={ 


Uhb 

^hb 


[-1 

► + < 

P 1 

► x < 

^ 1 

Vh 

> = < 

u + qz h - ry h j 
v + rx h - pz h 

Whb 


l ™ J 


i r J 


1 Zh ) 


w + py h ~ qx h J 


( 2 ) 


The translational acceleration at the rotor hub, in body axes, may be calculated in a similar manner. 
We have: 


where, 


V hb — V hb T ^ cgb ^ Vhb 


V hb = \ 


u + qz h - ry h 
v + rx h - pz h 
w + pyn - qx h 


(3) 


(4) 


Performing the cross product and summing, the acceleration vector at the rotor hub (body axes) is 
found to be: 
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( 5 ) 


{ u + qw - rv - ( q 2 + r 2 )xh + (pq ~ r)y h + (pr -f q)zh 
v + ru-pw + (pq + f)x h - (r 2 + p 2 )y h + (qr - p)z h > 
w + pv-qu+(pr - q)x h + (qr + p)Vh ~ (. V 2 + Q 2 ) z h 


To simplify the inclusion of the gravity terms in the analysis, the acceleration of gravity may be added 
to the hub acceleration vector at this point. This is similar to the way GenHel (ref. 7) deals with the 
gravity terms. The acceleration of gravity in body axes is found from: 


“ 1 

0 

0 


cos# 

0 

— sin# 

[ 0 ] 



—g sin# j 

0 

cos <p 

sin (p 


0 

1 

0 

0 

> = < 

9 

sin 4> cos # 

0 

— sin <p 

cos 4> 


sin# 

0 

cos# 

l 9 \ 


9 

cos 4> cos# J 


( 6 ) 


where # and (p are pitch and roll attitudes of the helicopter, respectively. Since gravity is an external 
force, the negative of this acceleration should be added to the hub acceleration which will then be used 
for calculating the inertial forces and moments. The augmented hub acceleration vector is: 



\ a hbx ] 

j it + qw - rv - ( q 2 + r 2 )xh + (pq - f)yh + ( pr - q)z h + g sin # 

Uhb = < 

1 1 
| a hby | 

> = < v + ru — pw + (pq + r)xh — ( r 2 + p 2 )yh + (<l r ~ fi) z h ~ 9 sin 0 cos # > 


( Q'hbz J 

[ w + pv - qu + (pr - q)xh + (qr + p)yh — ( p 2 + q 2 ) z h ~ 9 cos <p cos # 


To transform the velocity and acceleration vectors from the body axes to a fixed-shaft axes (fig. 2(b)), 
assuming that the shaft is tilted longitudinally through an angle ig (tilt back positive) and laterally 
through an angle (tilt right positive), the acceleration vector in body axes should be multiplied by 
the transformation matrix: 


[■ Ttut ] — 


COS ig 

sin i<p sin*# 
cos sin ig 


0 — sin ig 

cos {(f, sin cos ig 
— sin i^ cos i^ cos ig 


Therefore, the hub velocity and acceleration vectors in a fixed-shaft axes are: 



\ 

Uhfs 


^ hb 

Vhfs — < 

v hfs 

> — [Ttut] < 

Vhb 


Whfs 


k Whb 

I 

a hfsx ^ 


f 

a hbx 

O'hfs \ 

a hfsy | 

> — [Ttut] * 

O'hby 

1 

O’hfsz J 


K a hbz 


( Pfs 


\ V 

^ hfs ~ 

= ) 

> = [ Ttut ] ' 

* 


[ r f s . 


i r 


( 8 ) 


0 ) 


( 10 ) 


( 11 ) 
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Velocity and acceleration vectors at the blade element- In order to derive the inertial 
portion of the coupled flap-lag equations of motion using a Newtonian approach, the acceleration 
vector at an arbitrary blade element needs to be determined. The acceleration at the rotor hub in 
a fixed-shaft-axes system was derived previously. This acceleration vector has to be transformed to 
a rotating-shaft system and summed with the local acceleration vector. Referring to figure 2 (c), the 
hub acceleration in a rotating-shaft system may be calculated by rotating the fixed-shaft-frame- vector 
through an angle (7t/2 — ip) in the positive ^-direction. This is equivalent to multiplying by the 
rotation matrix: 


[Tfsrs] 


cos(7r/2 — ip) sin(7r/2 — ip) 0 
— sin(7r/2 — ip) cos(7r/2 — ip) 0 
0 0 1 


Simplifying the matrix and multiplying, we get: 


(12) 



sin^ 

cos ip 

0 ' 

f O'hfsx 

O'hrs — 

— cos ip 

sin ip 

0 

> a hfsy > 


0 

0 

1 

&hfsz 


{ dhfsx sin Ip + CLhfsy cos ip 1 

-dhfsx cos Ip + dhfsy sin ip > (13) 

Qhfsz J 


To calculate the contribution of the local flapping and lead-lag motion we first need to write the position 
vector of an arbitrary blade element in the rotating-shaft frame. Referring again to figure 2(a), the 
position vector in frame 1 is (note, that lead is positive): 


in frame 2 : 


in frame 3: 


r$ sin <5 
< r ( 5 cos 6 > 
0 


{ rg sin <5 
rg cos 8 + Ae 
0 


' 1 

0 

0 


| r$ sin 8 j I 

f rg sin 8 | 

0 

cos f3 

sin/5 

< 

rg cos 8 + Ae > = < 

cos /5 (rg cos 8 + Ae) ) 

0 

— sin/5 

cos /5 


{ o j 1 

1 — sin /5 (rg cos 8 + Ae) J 


and finally in the rotating-shaft frame: 


(14) 


T rs 


rg sin 8 

< cos f3 (rg cos 8 + Ae) + ep 
— sin (3 (rg cos 8 + Ae) 


(15) 
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The blade-element velocity vector in the rotating-shaft frame may then be calculated from: 


where the hub velocity vector in a rotating-shaft frame is: 

Uf s sin ip + v f s cos ip 
—Uf s cos ip T Vf s sin ip > 
w fs 

and the local derivative of the position vector is: 



V r . 


= { V r 


rsy 

V rS z 


— T rs — Vh rs + Vrs T 


X r r 


( 16 ) 


(17) 



rs 6 cos 8 

—f3 sin (3 (rs cos 6 + Ae) — rs 6 cos (3 sin S > 
—p3 cos (3 (rs cos 6 + Ae) + r 6 S sin (3 sin 6 


(18) 


The rotational velocity of the rotating-shaft axis is found by summing the rotations due to body and 
rotor motions, as follows: 


0 


sin ip 

cos ip 

0 ' 


' Pfs ' 


Pfs sin 'ip + qfs cos ip ) 

0 

> + 

— cos ip 

sin 'ip 

0 

< 

Qfs 

> = < 

~Pfs cos ip -P qf s sin 'ip 

l J 


0 

0 

1 


^ J 


r fs -n J 


MACSYMA was used for the derivation of the terms. Parameters in the MACSYMA derivations 
were given symbolic names representative of the notation used up to this point. For example, the 
blade-element velocity vector in the rotating-shaft frame (eq. 16) was developed as: 

VRS - RRSDI = VHRS + RRSDL + CROSS3(OMEGARS , RRS) 

where, 

RRSDI ~ r rs (r, rotating-shaft, derivative, inertial) 

RRSDL = f rs (r, rotating-shaft, derivative, local) 

and “CROSS3(X,Y)” is a simple MACSYMA macro that calculates the cross product of two 3-D vectors 
X and Y (included in app. B). The MACSYMA results are shown in appendix A. This velocity vector 
is later used, along with inflow and wind velocities, to calculate the aerodynamic forces and moments. 


The blade-element acceleration vector in a rotating-shaft frame: 


CLrs 


Q>rsx 

&rsy 

CLrsz 


— y r8 — ahrs + r rs + 2c u rs x f rs + tu rs x r rs + uj rs x (cu TS x r rs ) 


is found in a similar manner: 


( 20 ) 
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ARS = 


AHRS + RRSDDL + 2 * CROSS3{OMEGARS, RRSDL ) + CROSS?>{OMEGARSD , RRS) 

+CROS£3(OM£GARS, CROSS3(OMRGARS, RRS)) 
where, 

RRSDDL = r rs (r, rotating shaft, derivative, derivative, local) 

OMEGARSD = u> rs (omega, rotating shaft, derivative) 

MACSYMA results are again shown in appendix A. 

Inertial Terms 

Now that we have the inertial acceleration vector in the rotating-shaft frame, we can derive the 
inertial forces and moments and the inertial portion of the flap-lag equations of motion. 

Blade inertial shears at the lead-lag hinge, rotating-shaft frame— The elemental inertial 
force may be written, in rotating-shaft frame, as: 


A Fi rs ~ — Am a rs (21) 

The inertial force itself may then be written, in component form, as the sum of the elemental forces 
over the entire blade: 


F ' 

1 irsx 

j f £ Am &rsx 

F 

1 irsy 

/ — \ £ Am (Irsy 

L irsz j 

1 [ £ “Am a rsz 


( 22 ) 


where the summation is over all the blade-elements. Allowing MACSYMA to do the integration by 
letting: 


M — Y Am : Mass of one blade outboard of outer hinge 

FM = £ Am r$ : First mass moment of one blade about outer hinge 

we find expressions for the components of the inertial shear at the lead-lag hinge in rotating-shaft axes. 
The resulting equations are quite long and are, therefore, not included. In fact, from here on only the 
MACSYMA output for the flapping and lead-lag equations have been included. All the intermedaite 
equations, however, can be genetated by following the derivations provided. 

Blade inertial hub forces, rotating-shaft frame— The inertial forces at the hub in a rotating- 
shaft frame, for one blade, are the same as the inertial shears at the hinges. These were calculated 
before as: 


F%rsx 


£ Am drsx 1 

F 

j ^irsy 

► = < 

£ Am drsy j 

Dirsz 


£ “Am d rsz J 


(23) 


Blade inertial hub forces, fixed-shaft frame- The inertial hub forces in the fixed-shaft frame 
may be written, for one blade, as: 
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sin ip — cos ip 0 


Firsx 


F'irsx sin "Ip Fi rS y COS Ip 1 

cos ip sin ip 0 

0 0 1 


F 

r irsy 

Fjrsz 

> = < 

Firsx cos ip T F{ rS y sin ip , 
^ Firsz j 


( 24 ) 


Blade inertial moment at the lead-lag hinge, frame 1- In order to derive the inertial 
moment at the lead-lag hinge, in frame-1, the acceleration vector needs to be defined in frame 1 (refer 
to fig. 2(a)). Since the acceleration vectors in frames 1 and 2 are identical and frame 2 is simply a 
single rotation through the flapping angle (3 up from the rotating-shaft frame, we have: 




f &2x 


I a 2y 

> — 

{ a 2z 



1 0 0 

0 cos (3 — sin (3 
0 sin / 3 cos / 3 


O-rsx 


J &rsx | 


j" 

1 

a rsy j 


a rsy cos f3 - a rsz sin / 3 

l = < 

i 

drsz j 


1 a rsy sin / 3 + a rsz cos (3 J 


( &1 z 


Now, the moment at the lead-lag hinge may be written as: 

Mu l = yVjs] x (-Am Si) 


a l 

(25) 

(26) 


where, 



rg sin <5 
rg cos 6 
0 


> 


Therefore, the inertial moment at the lead-lag hinge is: 


(27) 


Mm — < 

Mu\x 


X^—Amr^aiz cos <5 ] 

Milly 

> = < 

E Am rg a iz sin 8 


K Muu , 


E~ Am rg (a\ y sin 6 — a\ x cos 8) J 


(28) 


Again allowing MACSYMA to do the integration, by additionally letting: 

SM = E A m r| : Second mass moment of one blade about the outer hinge 
we get expressions for the components of the inertial moment at the lead-lag hinge. 

The lead-lag equation of motion is simply the expression of the moment equilibrium about the 
lead-lag hinge ( 21 -axis) and will therefore only involve the Mm z component of the inertial moment 
at the lead-lag hinge. Also, the lead-lag hinge does not support any moment about the 21 -axis and 
therefore only the x\ and y\ components of the inertial moment are transferred down (except through 
the lead-lag dampers, which are considered independently). 

Blade inertial moments at the flapping hinge, frame 2— Inertial moments at the flapping 
hinge come from two sources, 1) the moment of the inertial shear forces at the lead-lag hinge, 2) the 
moment transferred through the lead-lag hinge. 

The inertial moment due to inertial shears at the lead-lag hinge may be written, in frame 2, as: 
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o ) ( F iflx | C A eF ifu ) 

M i/2 i = < Ae ^ x { Fi f iy l ^ 0 ) (29) 

0 J l Fifu j l -Ae F inx J 

where, 

' F iflx | r 1 0 0 1 C F irsx ' 

< Fifiy > = 0 cos /? -sin/? < F irsy > (30) 

Fifi z J _ 0 sin/3 cos (3 F irsz 

The inertial moment transferred through the lead-lag hinge may be written, in frame 2, as: 

Mmx 1 

Mi/22 — < Muiy > (31) 

0 J 

Then, the total inertial moment at the lead-lag hinge may be calculated as the sum of the two parts 
above, 

Mif 2 x ) Mmx + Ae Fifi z 1 

Mif2 = Mi f 21 + Mif22 — i Mif2y > = < Miliy > (32) 

Mif2z J k Fifi x J 

The flapping equation of motion is simply the expression of the moment equilibrium about the flapping 
hinge (x 2 -axis) and will therefore only involve the Mif 2x component of the inertial moment at the 
flapping hinge. Also, the flapping hinge does not support any moment about the .T 2 -axis and therefore 
only the y 2 and Z 2 components of the inertial moment are transferred down. 

Blade inertial hub moments, rotating-shaft frame- The inertial hub moment in the rotating- 
shaft frame, for one blade, is composed of two parts, 1) the moment of the inertial shear force at the 
flapping hinge, and 2) the moment transferred through the flapping hinge. 

The moments due to inertial shears at the flapping hinge may be written, in rotating-shaft frame, as: 

{ 0 Firsx 1 f Fi rsz 

Bp > X < Fi rS y / = \ 0 ► (33) 

0 y ^ Firsz J ^ ~€/3 Fi rsx 

The moments due to nonflapping moments at the flapping hinge may be written, in rotating-shaft 
frame, as: 

r 1 ° ° i r ° ] ( o 

M ihrS 2 = 0 cos/? sin/? < Mm y > = < M U i y cos/? - Ae F if i x sin/? > (34) 

_0 -sin/? cos/? J -Ae Fifi x j [ -Mm y sin/? - Ae F iflx cos (3 
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So the total inertial moment at the hub may be calculated as the sum of the two parts above, 


Mihrs — Mihrsl "F M \hrs2 — \ 


F irsz 

Muiy cos /? - Ae Fifi x sin/? 

M-iliy sin/? Ae cos/? e^g Fi rsx 


i 


( 35 ) 


Blade inertial hub moments, fixed-shaft frame- The total inertial moment at the hub, in 
the fixed-shaft frame, can now be calculated by rotating from rotating to fixed-shaft frame: 


Mi 


ihfs 


Multiplying, we get: 
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eR Fir cr cos ib + Mai,, cos/? 


(36) 


(37) 


Aerodynamic Terms 

The calculation of the aerodynamic forces and moments follows closely the implementation used 
by Howlett in GenHel (ref. 7). Specifically: 

1. Treatment of the blade segment aerodynamic force calculation is completely nonlinear. 

2. Bivariate maps, as a function of angle of attack and Mach number, are defined in the range -5 deg 
to +30 deg, allowing some accounting of blade stall effects. 

3. For angles of attack from +30 deg to +355 deg, values of section lift and drag coefficients are from 
C-81 equations. These, along with the low- angle maps, provide a complete coverage of angles of 
attack, allowing some treatment of the aerodynamic characteristics on the retreating blade side 
of the rotor disk. 

4. Simple sweep theory is used to modify the unyawed blade-element lift coefficients. 

5. Reynolds number effects, unsteady flow effects, and compressibility have been ignored. 

No attempt was made to improve the application of swept wing theory used to calculate elemental 
force components (tangential, radial, and normal) based on the lift and drag coefficients. The tables of 
coefficients themselves, however, and the lookup scheme used to access them were, of course, changed 
to represent the AH-64A Apache helicopter. Also, the different hinge sequence and the noncollocated 
hinges of the AH-64A required reformulation of the transfer of aerodynamic forces and moments. 
This in turn significantly modified the lead-lag and flapping equations as is detailed in the following 
derivations. For details of the lift and drag calculations using coefficients based on local angle of attack 
and Mach number, and the details of sweep theory used for enhancing the 2-D nature of the lift and 
drag coefficient data, refer to Howlett (ref. 7). 
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Blade aerodynamic shear forces, rotating-shaft frame- Let F ti$ , F ris , and F Pis be the 
tangential, radial, and perpendicular (normal) components of the elemental aerodynamic force at a 
given blade element (fig. 2(d)). Then the aerodynamic shear force may be written, in frame 1, as: 


cos# 

sin# 

0 ' 

[ -T.Ftu 

j f —Ft cos 8 F r sin# j 

— sin # 

cos# 

0 

ZF ri , 

> = < Ft sin # + F r cos # 
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0 

1 

{ ~EF Pi , 

1 1 j 


(38) 


The aerodynamic shear at the flapping hinge, at the lead-lag hinge, and at the hub are identical. 
Therefore, the hub aerodynamic shear in the rotating-shaft frame, for one blade, may be written as: 


or, 


1 0 0 
0 cos / 3 sin / 3 
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(39) 


F 1 

1 arsx 


—F t cos # + F r sin # ' 

< F 

' ± arsy , 

> — < 

Ft sin # cos P -j- F r cos # cos p — F p sin p > 

Farsz J 


[ —F t sin# sin/3 — F r cos# sin P — F p cos P 


(40) 


Blade aerodynamic shear forces, fixed-shaft frame- Rotating down to a fixed-shaft frame, 
we get: 
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(41) 


Blade aerodynamic moment about the lead-lag hinge, frame 1— The elemental aerody- 
namic force at each blade element, in frame- 1, is: 
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(42) 


the elemental aerodynamic moment at the lead-lag hinge, in frame-1, may then be written as: 



^ r 6 is s in # 1 

-F tls cos 6 + Fr is sin# ' 
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The blade aerodynamic moment about the lead-lag hinge may then be written as: 
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Mall = 

where the summation is over all the blade-elements. 

Blade aerodynamic moment at the flapping hinge, frame 2— The aerodynamic moment 
about the flapping hinge consists of two parts: 1) moment due to the shear force at the lead-lag hinge, 
and 2) moment due to the non-lead-lag component of moment about the lead-lag hinge. 

The moment due to aerodynamic shear force at the lead-lag hinge may be written, in frame 2, els: 

{ 0 1 ( —F t cos 8 + F r sin 8 ) i —A e F p 1 

Ae > x < Ft sin 8 + F r cos 6 > = < 0 > (45) 

0 J [ — F p j Ae Ft cos <5 — Ae F r sin 8 J 

The moment due to the non-lead-lag components of the moment about the lead-lag hinge may be 
written, in frame 2, as: 


z-r Si , F n, msS 

T,rs lt F pu sm 6 {. (44) 

£r&. F, 
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So, the total aerodynamic moment at the flapping hinge is: 


(46) 
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Blade aerodynamic hub moment, rotating-shaft frame- The aerodynamic hub moment 
consists of two parts, 1) the moment due to shear force at the flapping hinge, and 2) the moment due 
to the nonflapping moment at the flapping hinge. 

The moment due to shear force at the flapping hinge may be written as: 


or, 
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The moment due to the nonflapping moment about the flapping hinge may be written as: 


( 49 ) 
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or 
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So, the aerodynamic hub moment in a rotating-shaft frame, for one blade, is: 


(51) 
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Rotating to a fixed-shaft frame, we have: 
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Blade Restraint Terms 

Both lead-lag and flapping spring and damping constraints are dealt with as proportional terms. 
Referring to figure 2(e), the lead-lag spring and damper exert a moment of: 


Mm = I 


° . } 
k 6 s + k 6 s J 


(54) 


on the blade. The lead-lag spring and damper and the flapping spring and damper together exert a 
moment of: 


( Kp(3 + Kpp ] 

M r f2 = < 0 / (55) 

1 -K S 6- K h 8 J 

on the blade segment between the flapping and the lead-lag hinges (Ae in fig. 2(e)) in frame 2, and a 
moment of: 


15 



or 
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on the hub in the fixed-shaft axes. 

Lead-lag equation of motion— The lead-lag equation of motion is simply the solution of: 


Muiz T M a ii z T MfHz — 0 (58) 

for 8. This was done using MACSYMA with the “SIMPORG” macro (app. B) used for the factoring 
work. “SIMORG” basically takes the desired parts of the unsimplified equation from the initial stack 
and adds them to the final stack after performing algebraic and trigonometric simplifications (for more 
information refer to app. B). The resulting equation is given in appendix A. 

Flapping equation of motion— The flapping equation of motion is simply the solution of: 


f2x + M a f2 x + M r f2x — 0 . (59) 

for (3. This was also done using MACSYMA. The “SIMPORG” macro was again used for the factoring 
work, and the resulting equation is given in appendix A. 


Upgrades to Other Modules 

As in actual rotors, the total rotor outputs from the blade-element rotor module contain high 
frequency harmonics. Because of this harmonic variation, the original trim module supplied with 
FLYRT could not be used to trim the new model. The original trim module was designed to calculate 
the variation matrix only once and then use it throughout the trim procedure. A routine which would 
recalculate this matrix at every step was needed instead. Also, the original equations-of-motion module 
had to be modified to allow the averaging necessary to remove the harmonics of the forces and moments 
generated by the blade-element rotor module. 

A new trim module, based on the IMSL (ref. 13) subroutine ZXSSQ but following the general 
setup of the original FLYRT trim module, was therefore developed. Also, the equations-of-motion 
module was upgraded to provide the needed averaging during trim. The combination of the new trim 
and equations-of-motion routines were shown to be satisfactory for trimming BEMAP throughout the 
flight envelope. 

Further modifications of the equations-of-motion module were also necessary. In the map-type 
rotor module inertial effects are not taken into account explicitly. Instead, the polar moment of inertia 


16 



of the rotor is lumped with the other inertias in the equations-of-motion module. By contrast, the 
blade-element module explicitly accounts for the inertial forces and moments at the rotor. Therefore, 
the equations-of-motion module had to be modified to account for this difference. In effect, the rotor 
polar moment of inertia was replaced with the hub polar moment of inertia and the rotor inertias were 
accounted for explicitly at the rotor module itself. 


Apache Flight Data 

During the summer of 1990, an instrumented AH-64A Apache helicopter was used to perform a 
series of tests to collect data specifically for model validation purposes. The tests were conducted at 
the U.S. Army’s Airworthiness Qualification Test Directorate (AQTD) at Edwards Air Force Base in 
California, and were part of a broader effort (ref. 14) aimed at validating the new Handling Qualities 
Requirements for Military Rotorcraft (ADS-33C) (ref. 15). This database represents the best available 
source of flight-data for the AH-64A Apache. The only major shortcoming of the data is the fact that 
the helicopter used for the tests did not have an instrumented rotor and therefore rotor data were not 
collected. 

Both static and dynamic validation data were collected. Static (or trim) validation data were 
collected for various flight conditions from hover to 129 kts forward airspeed, several side- ward and 
rear-ward speeds, and a steady climb condition. In addition, some of the forward flight data were 
collected at two different aircraft C.G. positions to allow analysis of the effects of C.G. location on 
aircraft response. To collect trim data at each flight condition, the aircraft was flown to a stable trim 
at the desired flight condition and about 10 sec of data were collected at the trim state before moving 
on to the next run. 

The dynamic validation data consisted of doublets and frequency sweeps in all four control axes at 
hover, 60 and 120 kts forward airspeed. Data were collected with both DASE on and off. In order to 
reduce the effects of winds and gusts on the recorded responses, the flights were conducted in conditions 
of near zero winds. 


BEMAP Validation 

BEMAP was extensively validated by comparing its responses with the-flight data described above. 
The validation effort spanned both hover and forward flight and was conducted in both time and 
frequency domains. Some of the results of the validation work are presented below. 

Trim Validation 

The trim checks were conducted by setting up the model to match, as closely as possible, the 
configuration and flight condition of the available flight data. This included the use of flight-test values 
of aircraft weight, C.G. location, ambient temperature/density, and altitude. It’s usual practice to trim 
models at zero sideslip for low-speed flight and at zero-roll angle for high speed. Since the longitudinal 
and lateral components of the aircraft airspeed were available, however, it was decided to trim the 
model on roll angle at all speeds and treat sideslip angle as a fixed trim parameter. Therefore, four 
controls and two attitudes (pitch and roll) were used to trim the model in all cases. As discussed earlier 
a new trimming algorithm was developed to trim the new blade-element rotor. Using this routine, and 
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a simple averaging scheme in the equations-of- motion module (only during trim), it was possible to 
zero out the net forces and moments on the aircraft even with the harmonic nature of the rotor forces 
and moments. 

The checks were run for hover, 10, 20, 30, 40, 60, 80, 100, 120, and 129 kts and the results are 
shown in figures 3(a-g). As may be seen, the model trim values match flight data quite well. The 
differences are seen to be less than 10 percent in most cases. The largest veriations are in the 20-40 kts 
range which expands the transitional lift region where inflow modeling is less accurate. 

Figure 3(a) compares the predicted trim pitch attitude (both BEMAP and FLYRT) with flight 
data. The match between BEMAP and flight data is quite good close to hover and at high speed 
(>80 kts). Between 20 and 80 kts, however, the model first over and then under predicts the trim pitch 
attitude. The largest variation occurs at 40 kts where BEMAP underpredicts the actual pitch angle 
by about 2.9 deg. Everywhere else the variation is significantly smaller than this maximum. 

Figure 3(b) compares the predicted collective input required to trim with the flight values. BEMAP 
predicts the general trend of the collective variation with airspeed quite well. The model, however, 
consistently underpredicts the required collective by up to 10 percent. The variation is the largest at 
hover and high speed and is on average lower between 20 and 80 kts. The reason for this discrepancy at 
hover may be the simplified treatment of the rotor/ fusel age interaction incorporated into the fuselage 
model. The effect of the rotor down wash on the fuselage does not seem to be accounted for in a 
satisfactory manner, underpredicting the downforce on the fuselage. Note, that FLYRT does a slightly 
better job of predicting the collective input required to trim. 

Figure 3(c) compares the predicted longitudinal cyclic input to trim with the corresponding flight 
values. BEMAP is seen to duplicate the general trend in the migration of the longitudinal cyclic with 
airspeed quite well. The match is seen to be excellent below 20 kts (less than 1 percent difference) 
and quite good above 60 kts (around 3 percent or less). In the middle speeds (transitional lift region), 
however, the model overpredicts the longitudinal cyclic required to trim by up to 8 percent at one 
point. 

Figure 3(d) compares the predicted roll attitude with the actual flight values. The general trend 
of the roll angle variation with airspeed is duplicated very well. BEMAP slightly underpredicts the 
roll angle (by less than 1 deg) close to hover and then overpredicts it by a maximum of about 1 deg 
between 30 and 120 kts. 

Figure 3(e) compares the lateral cyclic required for trim with the flight values. Again, the general 
trend in the migration of the lateral cyclic with airspeed is duplicated quite well. Below 30 kts, BEMAP 
underpredicts the lateral cyclic required to trim by 1-2 percent. Above 30 kts, the model consistently 
overpredicts the lateral cyclic requirement by around 2 percent all the way to 120 kts. Overall, BEMAP 
results are seen to be closer to the flight values than FLYRT. 

Figure 3(f) compares the modeled pedal requirement for trim with flight data. It can be seen that 
the general variation with airspeed is not duplicated as well as it was for the other trim parameters. 
The difference between predicted and actual values, however, never exceeds 14 percent and is less than 
10 percent in most cases. The inadequate treatment of the effects of main- rotor downwash on the tail 
rotor is thought to be the cause of this discrepancy. 
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Finally, figure 3(g) shows the variation of the horizontal stabilator incidence angle with airspeed 
and shows that the model predicts the incidence angle quite well throughout. This indicates that 
an error in the horizontal stabilator angle is not responsible for the trim pitch attitude discrepancies 
between 20 and 80 kts. This however, does not mean that the horizontal stabilator is not the cause 
since the main rotor to horizontal tail inter ferance effect, and not the angle, may be the source of the 
errors. 


Time Domain Validation 

Similar to trim, the time-domain validations were conducted by setting up the model to match, as 
closely as possible, the configuration and flight conditions of the available flight data. This included the 
use of flight-test values of aircraft weight, C.G. location, ambient temperature/density, and altitude. 
No information, however, was available concerning moments and products of inertia of the test aircraft. 
Therefore, best estimates based on values given in the FLYRT documentation and related AH-64A work 
for an AH-64A helicopter configured similar to the test aircraft were used (table 1). 

Actual pilot-control inputs were used to drive the model in all cases. To avoid the initial abrupt 
control changes that would have been encountered due to mismatches in trim, the actual pilot controls 
were first converted to difference values. This simply consisted of subtracting the initial trim-control 
values from the control values at all subsequent times. These difference values were then added to the 
trim values of the controls calculated by the model. As a result, the lines depicting the model control 
inputs in all the plots presented in this section are parallel to the lines depicting the flight-test values. 
The bias between the two lines is an indication of how close the model trim values were to the actual 
flight-data values at the initiation of the data run. 

The flight data used for hover were collected with the DASE turned off whereas the flight data used 
for forward flight were collected with the DASE on. Since the validation of the airframe model was the 
only objective, control inputs were injected at the main- and tail-rotor swashplates for all maneuvers. 
This, in effect, bypassed the DASE module and prevented any DASE module errors from affecting the 
results of the comparisons. Note, however, that in the absence of the stabilizing effects of the DASE 
the model diverges much more rapidly. This point is explored later in this report. 

The test data were used as received from AEFA and no attempt was made to correct or remove 
instrument biases or recording errors. For each maneuver, control inputs in all four axes, angular 
rates, rotor rpm, body-axis velocities, normal acceleration, and horizontal stabilator incidence angle 
are shown and compared to flight data (in some cases engine torque is shown instead of stabilator 
incidence angle to help explain discrepancies). Note, that this large amount of data is provided for 
completeness and therefore not every plot is explicitly referred to in the following discussions. 

Hover 

Four maneuvers, a doublet in each control axis, were considered. The aircraft configuration data 
for the four maneuvers are presented in table 2. All the maneuvers were characterized by the pilot 
input in the main control being much larger than his inputs in the other three controls. Nevertheless, 
all the inputs were fed into the model to insure as close a match to the actual maneuver as possible. 
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In most cases, there is some small control activity in all the axes prior to the start of the maneuver. 
These were most likely made by the pilot to maintain or tighten the trim condition as the helicopter 
continually tends to depart from trim. Also, in some cases the intial rates at the beginning of the data 
are not zero and the pilot may have attempted to correct them using appropriate inputs prior to the 
actual initiation of the maneuver. Since the model rates are hard-zeros at trim, such corrective inputs 
can cause the model to diverge from the flight data even before the maneuver actually begins. This 
effect is especially pronounced at hover. To overcome this problem to some extent, a provision was 
added to allow the introduction of inputs to be delayed until a specific time after the start of the data 
(start of data is time — 0.0). The criterion for the selection of the modified start time was set to be the 
occurrence of zero rates (or as close to it as possible) in the flight data. In other words, the flight data 
for each maneuver were scanned and a time (prior to the beginning of the maneuver) was selected when 
all the angular rates were jointly as close to zero as possible. This time was then used as the start of 
pilot-input control of the model, the inputs being left at the trim values until then. As a consequence 
of this modification, the control input traces for the model are flat and do not follow the flight data 
between time = 0.0 and the modified start time. This modification was only used in hover. 

Finally, due to control freeplay in the actual aircraft, small inputs sometimes produced no flight-test 
response. Because there is no modeling of control freeplay, simulation responses may be expected to 
be faster with greater magnitude. A similar situation was noted by Ballin in his UH-60 modeling work 
(ref. 3). 

Lateral response— Figures 4(a-c) compare the responses of BEMAP with flight data for a left- 
right lateral doublet. The pilot inputs are depicted in figure 4(a). There is significant control activity 
prior to the start of the maneuver itself. These are especially noticeable in the lateral- cyclic and pedals. 
As discussed above, to prevent these extraneous inputs from affecting the match, the pilot-control-input 
injection into the model was delayed, in this case to time = 2.8 sec. Overall, BEMAP predicts the trim 
controls quite well in all axes except collective. As mentioned in the trim validation section, the main 
reason for the discrepancy in collective is thought to be insufficient modeling of the effects of rotor 
downwash impingement on the fuselage. Figure 4(b) shows that the roll-rate response of the aircraft 
is duplicated quite correctly by BEMAP. The yaw-rate response is also duplicated quite Well and the 
rotor rpm variation has the same character as the flight-test results. The pitch-rate response, however, 
initially goes the wrong way. The problem of predicting the pitch-to-roll and roll-to-pitch coupling 
responses will be discussed further, later in the report. Figure 4(c) shows that the components of 
body velocity are duplicated well initially. Finally, it seems that the incidence angle of the horizontal 
stabilator is overestimated by the model. However, according to the ADS- 10 document for the Apache 
(ref. 10), the horizontal stabilator should be at 25 deg in this flight condition. The reason for the lower 
flight data value is therefore not known. 

Longitudinal response- Figures 5(a-c) compare the responses of BEMAP with flight data for 
a forward-aft longitudinal doublet. Figure 5(a) shows that eventhough only longitudinal cyclic inputs 
were intended, there is significant activity in the lateral axes. Again, the model predicts the controls 
required to trim quite well except for collective. Note, that the injection of pilot-control into the model 
starts at time =1.0 sec. Figure 5(b) shows that the on- axis pitch rate is modeled quite well for the first 
part of the doublet. The match deteriorates, however, for the second part of the doublet. The off-axis 
roll response rate initially goes the wrong way in much the same way as the off-axis pitch response 
did in the lateral doublet case, suggesting a symmetry in the off-axis response error in hover. From 
figure 5(c) it may be seen that BEMAP predicts the normal acceleration well. Also, the variations in 
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the body lateral and longitudinal velocities are duplicated reasonably well initially. However, variations 
of the modeled horizontal tail incidence angle do not follow the flight data. Note, that the flight values 
of the horizontal tail incidence angle seem to exceed the 25 deg limit, implying a possible measurement 
error. 

Directional response- Comparisons of the responses of BEMAP with flight data for a left-right 
directional doublet are shown in figures 6(a-c). As may be seen from figure 6(a), pilot-input injection 
into the model does not occur until time = 2 sec. Figure 6(b) shows that yaw-rate response and 
rotor rpm variations are duplicated reasonably well. Also, BEMAP duplicates the aircraft roll rate 
reasonably well up to the second part of the doublet (9 sec). Finally, the pitch-rate response seems 
to follow the general trends of the actual response, with BEMAP and FLYRT responses being very 
similar. Figure 6(c) shows that the lateral translational velocity variation predicted by the model is 
significantly larger then the actual variation. The reason for this is probably the higher positive roll 
rate predicted by the model between 4 and 7 sec and the resulting higher roll attitudes. 

Heave response— Figures 7(a-c) compare BEMAP responses with flight data for an up-down 
collective doublet. Figure 7(a) depicts the control inputs for the maneuver. As may be seen from 
figure 7(b), all the model angular rate characteristics are similar to the flight data for the first 8-9 sec 
of the maneuver. The modeled rotor rpm response, however, shows significant discrepancy. The main 
reason for this error is the engine torque response, as depicted in figure 7(c). As may be seen, the 
engine response, as triggered by the collective input, is so fast that the model rotor not only does 
not droop, it actually overspeeds. No attempt, however, was made to modify the 701 engine module 
provided with FLYRT. Finally, figure 7(c) also shows that the model duplicates the normal acceleration 
response reasonably well, indicating proper implementation of the Pitt/Peters inflow model. 

60 lets 

Four maneuvers, a doublet in each control axis, were considered. The aircraft configuration data for 
the four maneuvers are presented in table 3. Again, all the maneuvers were characterized by the pilot 
input in the main control being much larger than his inputs in the other three controls. Nevertheless, 
all the inputs were fed into the model to insure as close a match to the actual maneuver as possible. 

Lateral response— Comparisons of BEMAP responses with flight data for a right-left lateral 
doublet at 60 kts are presented in figures 8(a-c). The pilot control inputs are shown in figure 8(a). As 
noted in the trim validation section, the pedal input required to trim is underpredicted by BEMAP. 
Identical results were obtained using FLYRT. The angular rate and rotor rpm responses are shown 
in figure 8(b). BEMAP predicts the roll-rate response very well for the first half of the doublet, but 
underpredicts the roll-rate magnitude during the second half. BEMAP also predicts the yaw-rate 
response quite well for the first 6 sec of the maneuver. The predicted pitch-rate response, however, 
is significantly in error. Note, that the discrepancy is quite similar to the discrepancy noted in hover 
for the lateral doublet maneuver. Finally, BEMAP predicts the general trend of rotor rpm variation 
quite well, though the magnitude of the response is underpredicted. Figure 8(c) compares BEMAP’s 
longitudinal and lateral components of body velocity, and normal acceleration with flight data and also 
shows the position of the horizontal stabilator throughout the maneuver. The longitudinal velocity 
compares well with the flight data. The predicted lateral-velocity response is somewhat different than 
flight, however, and the maximum lateral velocity encountered is slightly higher. Also, BEMAP does 
not predict the variation of the normal acceleration very well for this maneuver. Finally, the model 
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predicts the initial incidence angle of the horizontal stabilator within 1 deg. The match worsens with 
time, but it’s probably because the model, as a whole, is diverging from flight data. 

Longitudinal response— Figures 9(a-c) depict the validation results for a forward-aft longitudinal 
doublet at 60 kts. The control inputs are shown in figure 9(a). Figure 9(b) shows that the on-axis 
pitch-rate response is predicted quite well. Also, BEMAP seems to predict the general shape of the 
off-axis roll-rate response quite well. This is interesting since BEMAP’s prediction of the off-axis pitch- 
rate response to a roll doublet at 60 kts seemed to go the wrong way. Also, this is different than hover 
where the results showed a symmetry in the off-axes responses, by both being in error in a similar 
way. Figure 9(b) also shows that yaw coupling to pitch is not predicted well. Further, there is a 
discrepancy in the rotor rpm response even though the engine response is modeled fairly well, as shown 
in figure 9(c). Finally, figure 9(c) also shows that the general shape of the normal acceleration response 
is duplicated well. 

Directional response- Figures 10(a-c) compare the responses of BEMAP with flight-test results 
for a left-right directional doublet. Figure 10(a) shows the control inputs for the maneuver. Figure 10(b) 
indicates that the on axes yaw-rate response is modeled quite well throughout the maneuver. Also, the 
roll rate is matched quite well for the majority of the time. The pitch rate, however, is not predicted 
well. Figure 10(c) shows that the longitudinal component of body velocity remains fairly flat and this 
characteristic is modeled well. The variation of the lateral component of body velocity is also modeled 
well in magnitude, but there seems to be a significant delay in its actual buildup which is not duplicated 
well. Finally, the error seen in the normal acceleration is likely caused by the error in the pitch rate. 

Heave response— Figures ll(a-c) present comparisons of BEMAP responses with flight data for 
an up-down collective doublet. Figure 11(a) shows the control inputs for the maneuver. Figure 11(b) 
shows that all angular rates are predicted reasonably well. The rotor rpm variation, however, shows 
significant discrepancy. A look at figure 11(c) shows that the reason for this discrepancy, as in the 
hover collective doublet case, is the engine response. As may be seen, the engine-torque response as 
triggered by the collective input is so fast that the modeled rotor not only does not droop, it actually 
overspeeds. Finally, figure 11(c) shows that the normal acceleration is duplicated quite well, indicating 
that the dynamic inflow model seems to be valid. 

120 kts 

In general, BEMAP responses at 120 kts are not as good as those in hover and 60 kts. This is likely 
because at 120 kts the aerodynamic theory used for the rotor-blade-segment aerodynamics begins to 
be inadequate. Effects such as reverse flow and compressibility become more pronounced at these high 
speeds. Since such effects are not accounted for rigorously (or not at all in the case of compressibility) 
in the model, the model’s ability to predict the aircraft responses deteriorates. Again, four maneuvers, 
a doublet in each control axes, were considered. The aircraft configuration data for the four maneuvers 
are presented in table 4. 

Lateral response— Figures 12(a-c) show the results of comparing BEMAP responses with flight 
data at 120 kts for a right-left lateral doublet. Figure 12(a) shows that all the controls required to trim 
are predicted reasonably well except for collective. Interestingly, the discrepancy in collective seems 
quite similar to what was observed in the hover cases. In the 120 kts case, however, the rotor downwash 
impingement on the fuselage cannot be the cause of the discrepancy because the skew angle is quite 
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large. Figure 12(b) shows that the initial roll-rate response is duplicated quite well by the model, but 
it underpredicts the magnitude of the response to the second part of the doublet. Also, as in the hover 
cases, the pitch-rate response to a lateral input is not modeled well, and initially goes in the opposite 
direction. Figure 12(c) shows that neither BEMAP nor FLYRT seem to be very accurate in predicting 
the lateral and longitudinal components of body velocity, or the variation of the normal acceleration. 

Longitudinal response- Figures 13(a-c) compare the responses of BEMAP with flight data for a 
forward-aft longitudinal doublet at 120 kts. The control inputs are shown in figure 13(a). Figure 13(b) 
shows the on-axis pitch response to be duplicated quite well. More interestingly, the initial off-axis 
roll-rate response of the model seems to follow the general trend of the actual response. Specifically, 
the response does not seem to go in the opposite direction, as it does in the case of the off-axis pitch 
response to a roll input. Since the same kind of results were noted at 60 kts, this further suggests that 
there is some asymmetry in the ability to duplicate the two off- axes responses in forward flight. This 
asymmetry might provide some clues as to the cause of the discrepancies and why none of the current 
models are capable of duplicating the off-axis response (as discussed in ref. 2). 

Directional response— Figures 14(a-c) show BEMAP responses to a directional doublet at 
120 kts. The pilot controls are shown in figure 14(a). Figure 14(b) indicates that the BEMAP yaw-rate 
response duplicates flight quite well. BEMAP pitch and roll rates, however, show significant discrep- 
ancy. Note, that FLYRT does a better job matching the flight-roll rate. Finally, figure 14(c) shows that 
both models fall short of matching either the velocity components or the normal acceleration observed 
in flight. 

Heave response— Figures 15(a-c) depict BEMAP responses to a collective doublet at 120 kts. 
Pilot inputs are shown in figure 15(a). Figure 15(b) shows that BEMAP duplicates the general trends of 
the angular-rate responses of the aircraft. Again, however, the fast engine-module response to collective 
causes the model rpm responses to be incorrect. Finally, as figure 15(c) shows, BEMAP fails to match 
the lateral component of velocity observed in flight. 


Effects of Bypassing the Date Model 

As mentioned earlier, in the absence of the stabilizing effects of the DASE the model diverges 
much more rapidly. Introducing pilot inputs at the swashplate, as was done here to concentrate on the 
airframe model validation, in effect removes the stabilizing effects of the modeled DASE. Moreover, since 
the test aircraft’s DASE could have been responding to conditions not duplicated by the simulation, 
significant discrepancies could result. To demonstrate this effect, the lateral maneuver at 120 kts 
discussed previously with pilot inputs introduced at the swashplate was repeated with the modeled 
DASE turned on and pilot inputs introduced at the stick. The results are shown in figure 16. As 
may be seen, the new results show significant improvement over the original run in the on-axis. The 
roll-rate response is seen to match the flight data much better, especially for the second half of the 
doublet. As expected, the effect of feedback is to mask errors at low frequencies. Therefore, leaving the 
DASE model on and introducing pilot inputs at the stick for all the runs would have led to much closer 
results. However, it would have masked errors in the airframe modeling which need to be identified 
and addressed. 
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Frequency Domain Validation 


Flight-data-based frequency responses for hover were already available from Schroeder et. al. 
(ref. 1). Frequency responses for forward flight, however, were not available and had to be generated as 
part of the BEMAP validation effort. As in reference 1, CIFER^ was used along with frequency sweep 
data collected at AEFA (discussed earlier) to generate the frequency responses. Briefly, the process 
was as follows: 

1. All available forward-flight-sweep data were processed through FRESPID and MISOSA using a 
single representative window size. FRESPID is a CIFER^ subprogram which uses the chirp- 
z transform (an advanced fast-fourier transform algorithm) and overlapping hanning windows 
(used to provide low spectral variance through averaging) to calculate magnitude, phase, and 
coherence data from frequency sweep time histories. MISOSA is another CIFER^ subprogram 
which processes FRESPID results to remove the corrupting effects of correlated inputs in the 
off-axis controls. In effect, the output of MISOSA consists of “conditioned” frequency responses 
that represent the relationship of the output to the input with the contribution of all other inputs 
removed. 


This preliminary processing was performed to insure that the highest possible coherence could 
be achieved in the final-frequency response pairs. 

2. For each axis, the two/three best runs were selected (based on the coherence plots obtained from 
MISOSA in step 1), concatenated, and reprocessed through FRESPID using five windows. The 
windows used were 10, 20, 30, 35, and 40 sec long and were sized to provide the best combination 
of low-frequency coverage and maximized spectral averaging. 

3. The results of step 2 were then processed through MISOSA to eliminate the contaminating effects 
of correlated off-axis inputs. 

4. Finally, the results of step 3 were processed through COMPOSITE, which is the CIFER H sub- 
program used to optimize the data obtained using various window sizes. 


Only DASE-off data were used for all the frequency-domain comparisons. In the absence of stability 
augmentation the pilot often has to use significant off-axis control inputs in order to maintain safe 
flying attitudes during the sweep. Consequently, the coherence of the frequency response pairs are 
poor in some cases, especially in the off-axis responses. 


Ideally, it would have been preferable to drive the model with the same frequency-sweep inputs 
and process the model responses through CIFER^ to generate frequency response pairs (magnitude 
and phase plots). However, since the open-loop model is unstable and would diverge long before the 
90 sec duration of a typical sweep, running frequency sweeps through the model would have required 
additional effort to introduce additional stability loops as in reference 3. Instead, a simple numerical 
linearization technique was used to generate linear state-space models which could then be compared 
with the flight data. The stability derivative generation routine delivered with FLYRT was modified for 
this task. In particular, it was modified to perform a two-sided linearization to provide more accurate 
results, i.e., each derivative was calculated as: 



y(x Q + Ax) - y(x Q - Ax) 
2Ax 


X=Xo 
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Hover 


Preliminary hover results were briefly documented in reference 16 and are expanded upon here. As 
in reference 16, the identified time delays from reference 1 were used with BEMAP 6-DOF results to 
partially compensate for the limitations of 6-DOF analysis. Table 5(a) provides a comparison of the 
stability derivatives identified by Schroeder et. al. (ref. 1) with the derivatives obtained from BEMAP. 
As may be seen, there is general agreement between the model and the identified on-axis rotational- 
stability derivatives. The pitch to roll (Af p ) and roll to pitch ( L q ) off-axis derivatives, however, do not 
agree with the flight-identified values. In fact, the BEMAP-based values for both M p and L q are of the 
same magnitude but opposite sign compared to the identified values. This is not surprising considering 
the time domain comparisons discussed earlier. 

Moving to control derivatives, examination of the results in table 5(b) shows good agreement for the 
on-axis rotational pitch (Mj on ), roll ( Li at ), and yaw ( N pe d ) control derivatives. The on-axis translational 
longitudinal ( Xion ), lateral (Yi at ), and heave (Z co i) control derivatives also show general agreement with 
the flight values. 

The closeness of the match between BEMAP responses and flight data can be better seen by compar- 
ing the actual frequency response pairs instead of derivative values! These are shown in figures 17(a-d). 
In figure 17(a), BEMAP roll-rate and lateral- velocity responses are compared with flight data. As can 
be seen, BEMAP responses match the flight data quite well, especially in the frequency range between 
0.8 and 7 rad/sec. Furthermore, BEMAP responses show slight improvement over FLYRT for both 
cases. Figure 17(b) compares BEMAP’s pitch-rate and longitudinal-velocity responses with flight data. 
Again, the figure shows that the responses are quite close. The coherence of the pitch-rate response 
results is quite good and the match is seen to be good in a wider frequency range (0.5-10 rad/sec) 
compared to roll. However, the unsatisfactory coherence of the longitudinal velocity results at very low 
frequencies (less than 0.5 rad/sec) makes them unreliable. Figure 17(c) compares BEMAP yaw rate 
and heave velocity responses with flight data. The responses again compare favorably and BEMAP and 
FLYRT seem to be identical in high frequencies (this should be expected since BEMAP uses the same 
tail- rotor module). At low frequencies, however, FLYRT seems to do better in predicting the yaw-rate 
response. Finally, figure 17(d) depicts the poor match between BEMAP off-axes responses and flight 
data and shows that FLYRT responses are the same way. Note, that the match in magnitude is signif- 
icantly better than the match in phase, which reflects the time domain effect of magnitude matching 
but with wrong phasing (opposite), and the derivatives M p and L q having the correct magnitude but 
the wrong sign. 


60 kts 

Tables 5(a-b) provide the stability and control derivatives obtained from BEMAP and FLYRT at 
60 kts. Since state-space models based on flight data were not identified in forward flight, the model 
derivatives can not be evaluated and are only provided for completeness. The quality of the BEMAP 
responses can, however, be evaluated by comparing Bode plots of the responses of a BEMAP-based 
6-DOF linear model (the same identified time delays used in hover were used here due to a lack of more 
appropriate values) with frequency response pairs obtained from flight data. Figures 18(a-c) show 
such a comparison. Figure 18(a) shows the roll- and pitch-rate responses to lateral and longitudinal 
inputs respectively. As may be seen, BEMAP does quite well in predicting the roll- rate response 
and significantly improves on FLYRT. The two models both match the pitch-rate data reasonably 
well. Figure 18(b) shows that BEMAP provides noticeable improvement over FLYRT in predicting 
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the yaw-rate response to pedals. Both models, however, have difficulty predicting the vertical- velocity 
response to a collective input. Finally, figure 18(c) corroborates the results obtained in the time domain 
regarding the asymmetry of the off-axis prediction in forward flight. As may be seen, BEMAP does 
fairly well in predicting the off-axis roll-to-pitch response at frequencies above 0.7 rad/sec, showing 
noticeable improvement over FLYRT. It still, however, does poorly in predicting the pitch response to 
a lateral input, as does FLYRT. 


120 kts 

The frequency response pairs generated from the 120 kts flight data had the lowest overall coherence 
of all the data. As a result only the on-axes responses had sufficient coherence to provide reasonable 
reliability. Therefore, only on-axes responses are presented here. Table 6(a-b) present the 6-DOF- 
stability derivatives from BEMAP and FLYRT. As in the 60 kts case, flight-identified state-space 
models were not available, so the derivatives can not be compared to flight data. They are, nevertheless, 
presented for completeness. Also, the hover flight-identified time delays were again used due to a lack 
of more appropriate values. Figure 19(a) depicts the on-axes lateral and longitudinal responses. As 
may be seen, BEMAP does quite well in both cases and shows improvement over FLYRT. Figure 19(b) 
depicts the on-axes directional and heave responses. Again, BEMAP responses are seen to match the 
flight data reasonably well and to show improvement over FLYRT. 


Concluding Remarks 

A blade-element mathematical model for the AH-64A Apache Advanced Attack Helicopter 
(BEMAP) has been developed by the Aerofiightdynamics Directorate of the U.S. Army ATCOM. 
BEMAP is based on the MDHS FLYRT model of the AH-64A, but incorporates a blade-element ap- 
proach in its main-rotor module. This approach treats the aerodynamic and especially the inertial 
forces and moments more rigorously and removes the dependency on pregenerated maps inherent in 
the original FLYRT main-rotor module. 

Results of the BEMAP validation effort described in this report indicate that: 

1. Responses compare favorably with FLYRT, indicating that the blade-element rotor module has 
been derived and implemented correctly. 

2. On-axis responses match hover and forward-flight flight data well as indicated by the time and 
frequency domain comparisons. 

3. Off-axes results match flight data poorly at hover for both pitch-to-roll and roll-to-pitch responses. 
This may be seen in both time and frequency domains and by comparing the off-axes stability 
derivatives with those identified from flight. 

4. Off-axes correlation results for forward flight are mixed. The aircraft’s roll response to a pitch 
input is predicted fairly well whereas the pitch response to a roll input is not. This difference can 
be seen in both the time and the frequency domain results. 

Improvements in our AH-64A simulation capability, including a solution to the off-axes discrepancies, 
can now be attempted by including enhanced aerodynamics and explicit treatment of compressibility, 
blade stall, and reverse flow in the blade-element rotor module. 
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Appendix A 



Blade-element velocity vector in rotating-shaft frame 


>' rsr — 

COS(PSI) VHFS + SIN(PSI) UHFS + (COS(BETA) (COS(DELTA) RDELTA + DELE) 

+ EBETA) (OMG - RFS) - SIN(BETA) (SIN(PSI) QFS - PFS COS(PSI)) 

dDELTA 

(COS (DELTA) RDELTA + DELE) + COS (DELTA) RDELTA 

dT 

V — 

* rsy 

SIN(PSI) VHFS - COS(PSI) UHFS + SIN (DELTA) RDELTA (RFS - OMG) 

- SIN (BETA) (- COS(PSI) QFS - PFS SIN(PSI)) (COS(DELTA) RDELTA + DELE) 

dBETA 

- SIN (BETA) (COS (DELTA) RDELTA + DELE) 

dT 

dDELTA 

- COS (BETA) SIN (DELTA) RDELTA 

dT 

v — 

' rs: — 

WHFS + (COS(PSI) QFS + PFS SIN(PSI)) 

(COS(BETA) (COS(DELTA) RDELTA + DELE) + EBETA) 
dBETA 

- COS (BETA) (COS (DELTA) RDELTA + DELE) 

dT 

+ SIN (DELTA) (PFS COS(PSI) - SIN(PSI) QFS) RDELTA 

dDELTA 

+ SIN (BETA) SIN (DELTA) RDELTA 

dT 


Blade-element acceleration vector in rotating-shaft frame 

Q r sx — 

dOMG dRFS 

(COS (BETA) (COS(DELTA) RDELTA + DELE) + EBETA) ( ) 

dT dT 
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+ (OMG - RFS) (SIN (DELTA) RDELTA (RFS - OMG) 

- SIN(BETA) (- COS(PSI) QFS - PFS SIN(PSI)) (COS (DELTA) RDELTA + DELE)) 

dBETA 

+ 2 ((- SIN(BETA) (COS(DELTA) RDELTA + DELE) 

dT 

dDELTA 

- COS (BETA) SIN (DELTA) RDELTA) (OMG - RFS) 

dT 

dDELTA 

+ (SIN(PSI) QFS - PFS COS(PSI)) (SIN(BETA) SIN(DELTA) RDELTA 

dT 

dBETA 

- COS(BETA) (COS(DELTA) RDELTA + DELE))) 

dT 

+ (SIN(PSI) QFS - PFS COS(PSI)) ((COS(PSI) QFS + PFS SIN(PSI)) 

(COS (BETA) (COS (DELTA) RDELTA + DELE) + EBETA) 

+ SIN (DELTA) (PFS COS(PSI) - SIN(PSI) QFS) RDELTA) 

dQFS dPSI dPSI 

- SIN(BETA) (SIN(PSI) + COS(PSI) QFS + PFS SIN(PSI) 

dT dT dT 

2 

dPFS d DELTA 

COS(PSI)) (COS(DELTA) RDELTA + DELE) + COS(DELTA) RDELTA 

dT 2 

dT 

dDELTA 2 

- SIN(DELTA) (— ) RDELTA + AHFSX SIN(PSI) + AHFSY COS(PSI) 

dT 


&rsy — 

dRFS dOMG 

SIN (DELTA) RDELTA ( ) + 2 

dT dT 

dDELTA 

(COS(DELTA) RDELTA (RFS - OMG) + (- COS(PSI) QFS - PFS SIN(PSI)) 

dT 

dDELTA dBETA 
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(SIN (BETA) SIN (DELTA) 


RDELTA - COS (BETA) 


(COS(DELTA) RDELTA + DELE))) + ((COS (BETA) (COS (DELTA) RDELTA + DELE) + EBETA) 
(OMG - RFS) - SIN(BETA) (SIN(PSI) QFS - PFS COS(PSI)) 

(COS (DELTA) RDELTA + DELE)) (RFS - OMG) 

+ (- COS(PSI) QFS - PFS SIN(PSI)) ((COS(PSI) QFS + PFS SIN(PSI)) 

(COS(BETA) (COS(DELTA) RDELTA + DELE) + EBETA) 

+ SIN (DELTA) (PFS COS(PSI) - SIN(PSI) QFS) RDELTA) 

dQFS dPSI dPSI 

- SIN (BETA) (- COS(PSI) + SIN(PSI) QFS - PFS COS(PSI) 


dPFS 

SIN(PSI) ) (COS(DELTA) RDELTA + DELE) 

dT 

2 

d BETA 

- SIN(BETA) (COS (DELTA) RDELTA + DELE) 

2 

dT 

dBETA 2 

- COS (BETA) ( ) (COS (DELTA) RDELTA + DELE) 

dT 

2 

d DELTA dDELTA 2 

- COS (BETA) SIN (DELTA) RDELTA - COS(BETA) COS (DELTA) ( ) RDELTA 

2 dT 

dT 


dBETA dDELTA 

+ 2 SIN(BETA) SIN (DELTA) RDELTA + AHFSY SIN(PSI) - AHFSX COS(PSI) 

dT dT 


U j rsz — 

(COS(PSI) QFS + PFS SIN(PSI)) (SIN (DELTA) RDELTA (RFS - OMG) 

- SIN(BETA) (- COS(PSI) QFS - PFS SIN(PSI)) (COS(DELTA) RDELTA + DELE)) 
+ (PFS COS(PSI) - SIN(PSI) QFS) ((COS(BETA) (COS(DELTA) RDELTA + DELE) 
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+ EBETA) (OMG - RFS) - SIN(BETA) (SIN(PSI) QFS - PFS COS(PSI)) 
(COS (DELTA) RDELTA + DELE)) + 2 ((COS(PSI) QFS + PFS SIN(PSI)) 
dBETA 

(- SIN (BETA) (COS (DELTA) RDELTA + DELE) 

dT 


dDELTA dDELTA 

- COS (BETA) SIN (DELTA) RDELTA) + COS (DELTA) 

dT dT 

dQFS dPSI 

(PFS COS(PSI) - SIN(PSI) QFS) RDELTA) + (COS(PSI) SIN(PSI) QFS 

dT dT 


dPSI dPFS 

+ PFS COS(PSI) + SIN (PSI) ) (COS (BETA) (COS (DELTA) RDELTA + DELE) 

dT dT 


2 

d BETA 

+ EBETA) - COS(BETA) (COS (DELTA) RDELTA + DELE) 

2 

dT 


dBETA 2 

+ SIN (BETA) ( ) (COS (DELTA) RDELTA + DELE) 

dT 


dQFS dPSI dPSI 

+ SIN (DELTA) (- SIN (PSI ) COS(PSI) QFS - PFS SIN(PSI) 


dT 


dT 


dT 


2 

dPFS d DELTA 

+ COS(PSI)) RDELTA + SIN (BETA) SIN (DELTA) RDELTA 

dT 2 


dT 


dDELTA 2 

+ SIN (BETA) COS (DELTA) (- ) RDELTA 

dT 

dBETA dDELTA 

+ 2 COS (BETA) SIN (DELTA) RDELTA + AHFS2 

dT dT 
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Lead-lag equation of motion 

6 = 


( (COS (BETA)* ’DIFF(RFS ,T , 1 )+( 1-COS (BETA) ** 2) *C0S( DELTA) *SIN (DELTA) * 

1 RFS**2+( ( ( 2*SIN ( BETA )*C0S( DELTA )**2-SIN( BETA) )*C0S(PSI)-2*C0S(B 

2 ETA) *SIN (BETA) *C0S (DELTA) *SIN(DELTA) *SIN (PSI ) ) *QFS+ (2*SIN (BETA) 

3 *COS (DELTA) **2-SIN (BETA) ) *PFS*SIN (PSI )+2*C0S (BETA) *SIN (BETA) *C0 

4 S (DELTA )*SIN (DELTA) *PFS*CQS (PSI )+(2*CQS (BETA) **2-2 )*C0S (DELTA)* 

5 SIN (DELTA ) *0MG-2*SIN ( BETA ) * ' DIFF (BETA , T , 1 ) *C0S (DELTA) **2) *RFS+S 

6 IN(BETA)*SIN(PSI)*'DIFF(QFS,T,1)+((C0S(BETA)-2*C0S(BETA)*C0S(DE 

7 LTA) **2)*C0S (PSI) *SIN(PSI)+ (-COS (BETA) **2-1) *C0S (DELTA )*SIN (DEL 

8 TA) * COS ( PSI )**2+C0S( BETA) **2*COS (DELTA )*SIN (DELTA) )*QFS**2+(SIN 

9 (BETA) *C0S (PSI )*’ DIFF (PSI, T, 1 )+ ( (-2*C0S(BETA) **2-2) *COS (DELTA)* 

: SIN (DELTA) *PFS*C0S (PSI )+2*C0S( BETA) *SIN (BETA )*C0S (DELTA )*SIN(DE 

; LTA) *0MG+2*C0S (BETA) * »DIFF(BETA,T, l)*C0S(DELTA)**2)*SIN(PSI)+(4 

< *COS (BETA) *COS (DELTA )**2-2*C0S (BETA) )*PFS*C0S (PSI )**2+( (SIN (BET 

= A)-2*SIN (BETA) *COS (DELTA) **2)*0MG+2* 'DIFF (BETA , T, 1) *COS (DELTA) * 

> SIN (DELTA) ) *COS (PSI) + (COS (BETA) -2*C0S (BETA) *C0S (DELTA) **2)*PFS) 

? *QFS+SIN (BETA) *PFS*SIN (PSI)* 'DIFF (PSI ,T, l)+( (2*C0S(BETA) *C0S(DE 

« LTA) **2-C0S (BETA) ) *PFS**2*C0S(PSI)+ ( (SIN(BETA) -2*SIN (BETA)* COS ( 

1 DELTA) **2) *0MG+2* ' DIFF (BETA , T , 1 ) *C0S (DELTA ) *SIN (DELTA) ) *PFS ) *SI 

2 N (PSI) + (COS (BETA) **2+1 )*COS (DELTA) *SIN (DELTA) *PFS**2*C0S(PSI)** 

3 2+( ( -2*C0S ( BETA) *SIN( BETA )*COS (DELTA )*SIN (DELTA)*0MG-2*C0S(BETA 

4 )*'DIFF(BETA,T, 1)*C0S(DELTA)**2)*PFS-SIN(BETA)* ’DIFF(PFS,T, 1))* 

5 COS (PSI) -COS (DELTA)*SIN (DELTA)*PFS**2-C0S (BETA) * 'DIFF(0MG ,T , 1 ) + 

6 ( 1-COS (BETA)**2) *COS (DELTA) *SIN (DELTA) *0MG**2+2*SIN(BETA) * 'DIFF 

7 (BETA ,T,1)*C0S (DELTA) **2*0MG-('DIFF(BETA,T,1))**2*C0S(DELTA)*SI 

8 N (DELTA) ) *SM+ ( COS (DELTA) *EBETA+COS ( BETA) *DELE*COS (DELTA ) ) *FM* ' D 

9 IFF (RFS , T , 1 ) + (-COS (BETA) *SIN (DELTA) *EBETA-COS (BETA) **2*DELE*SIN 

: (DELTA) ) *FM*RFS**2+( ( (-SIN (BETA) *SIN (DELTA) *EBETA-2*C0S (BETA) *S 

; IN (BETA) *DELE*SIN (DELTA) ) *FM*SIN (PSI ) +SIN (BETA ) *DELE*COS (DELTA) 

< * FM*COS (PSI ) )*QFS+SIN( BETA) *DELE*COS (DELTA )*FM*PFS*SIN (PSI)+(SI 

= N (BETA ) *SIN (DELTA ) *EBETA+2*C0S (BETA) *SIN (BETA ) *DELE*SIN (DELTA ) ) 

> *FM*PFS*C0S (PSI) +(2*C0S (BETA )*SIN(DELTA)*EBETA+2*C0S (BETA )**2*D 
? ELE*SIN(DELTA))*FM*0MG-2*SIN(BETA)*'DIFF(BETA > T,1)*DELE*C0S(DEL 

0 TA)*FM)*RFS+ (SIN (BETA) *DELE*C0S (DELTA) *FM*SIN(PSI)-SIN (BETA) *SI 

1 N (DELTA ) *EBETA*FM*COS ( PSI ) ) * 'DIFF (QFS , T , 1 ) + ( (-CDS (DELTA) *EBETA- 

2 COS ( BETA) *DELE*C0S (DELTA) ) *FM*COS (PSI ) *SIN (PSI ) + (-COS (BETA) *SIN 

3 (DELTA) *EBETA-COS (BETA) **2*DELE*SIN (DELTA) )*FM*C0S (PSI ) **2+ (COS 

4 (BETA) **2-1 )*DELE*SIN (DELTA) *FM)*QFS**2+( (SIN (BETA) *SIN(DELTA)* 

5 EBETA*FM*SIN (PSI ) +SIN (BETA) *DELE*C0S (DELTA) *FM*C0S (PSI ) ) * 'DIFF( 

6 PSI , T , 1 ) + ( (-2*C0S (BETA) *SIN (DELTA) *EBETA-2*C0S (BETA) **2*DELE*SI 

7 N (DELTA ) ) *FM*PFS*COS (PSI) + ( SIN (BETA) *SIN (DELTA ) *EBETA+2*C0S ( BET 

8 A ) *SIN (BETA) *DELE*SIN (DELTA) ) *FM*0MG+2*C0S (BETA ) * ' DIFF (BETA , T , 1 

9 ) *DELE*COS (DELTA) *FM) *SIN (PSI) + (2*C0S (DELTA) *EBETA+2*C0S (BETA) * 

: DELE*C0S (DELTA) ) *FM*PFS*C0S (PSI) **2+ (2* 'DIFF (BETA , T , 1) *DELE*SIN 

; (DELTA) +FM-SIN (BETA) *DELE*C0S (DELTA) *FM*OMG) *COS(PSI )+ (-COS (DEL 

< TA) *EBETA-COS (BETA) *DELE*COS (DELTA) ) *FM*PFS) *QFS+ ( SIN(BETA) *DEL 

= E*COS (DELTA) *FM*PFS*SIN (PSI) -SIN (BETA) *SIN (DELTA) *EBETA*FM*PFS* 

> COS(PSI))*'DIFF(PSI,T,l)+( ( COS (DELTA ) *EBETA+C0S (BETA) *DELE*COS ( 


33 



? DELTA) ) *FM*PFS**2*C0S (PSI ) -SIN ( BETA )*SIN (DELTA) *EBETA*FM* 'DIFF( 
3 PFS , T , 1 ) + ( 2* ' DIFF ( BETA , T , 1 ) *DELE*SIN (DELTA ) +FM-SIN ( BETA ) ♦DELE* C 

1 OS (DELTA) *FM*OMG ) *PFS+ ( AHFSY*COS ( BETA) *SIK (DELTA) -AHFSX+COS (DEL 

2 TA) ) *FM) *SIN (PSI ) + (COS (BETA) *SIN (DELTA) *EBETA+COS (BETA) **2*DELE 

3 *SIN (DELTA) ) *FM*PFS**2*C0S (PSI ) **2+ (-SIN (BETA) *DELE*COS (DELTA) * 

4 FM* 'DIFF (PFS , T, l)+( ( -SIN ( BETA) *SIN( DELTA )*EBETA-2*C0S( BETA )*SIN 

5 ( BETA ) *DELE*SIN (DELTA ) ) *FM*OMG-2*COS ( BETA ) * * DIFF ( BETA , T , 1 ) *DELE 

6 *COS (DELTA ) *FM)*PFS+ ( -AHFSX*COS ( BETA) *SIN ( DELTA) -AHFSY*COS (DELT 

7 A) ) *FM) *COS (PSI ) + ( -COS (BETA) *SIN (DELTA ) *EBETA-DELE*SIN (DELTA ) ) * 

8 FM*PFS**2+ ( -COS (DELTA ) *EBETA-COS (BETA) *DELE*COS (DELTA) ) *FM* 'DIF 

9 F ( OMG , T , 1 ) + ( -COS (BETA ) *SIN (DELTA) *EBETA-COS(BETA) **2*DELE*SIN (D 

: ELTA ) ) *FM*0MG**2+2*SIN ( BETA) * ’ DIFF ( BETA , T , 1) *DELE*C0S (DELTA) *FM 

; *OMG-MRLlZ-MALlZ+(-('DIFF(BETA,T,l))**2*DELE-AHFSZ*SIN(BETA))*S 
< IN (DELTA)*FM)/SM 


Flapping equation of motion 

8 — 


( ( SIN ( BETA) *COS (DELTA )*SIN (DELTA)* 'DIFF(RFS , T, l) -COS (BETA) *SIN (BET 

1 A)*COS(DELTA) **2*RFS**2+( ( (2*C0S( BETA )**2-l)*COS (DELTA )**2*SIN ( 

2 PSI ) +C0S (BETA) *COS (DELTA) *SIN (DELTA) *C0S ( PSI ) ) +QFS+COS (BETA) *C0 

3 S (DELTA) *SIN (DELTA)*PFS*SIN (PSI )+( l-2*COS (BETA) **2)*C0S (DELTA)* 

4 *2*PFS*C0S (PSI) +2*C0S (BETA) *SIN(BETA)*C0S(DELTA)**2*0MG+2*SIN(B 

5 ETA) *C0S(DELTA)**2*’ DIFF (DELTA ,T, 1) )*RFS+(C0S(DELTA)**2*C0S(PSI 

6 )-C0S (BETA) *COS (DELTA) *SIN (DELTA) *SIN (PSI))* 'DIFF(QFS,T,l)+(COS 

7 (BETA) *SIN (BETA) *COS (DELTA ) **2*SIN (PSI) **2+SIN( BETA) *C0S (DELTA) 

8 *SIN(DELTA) *C0S(PSI)*SIN(PSI))*QFS**2+ ((-COS (DELTA) **2*SIN(PSI) 

9 -COS (BETA) *COS( DELTA ) *S IN (DELTA) * COS (PSI) ) * 'DIFF (PSI , T , 1 ) +SIN(B 

: ETA) *COS (DELTA) *SIN(DELTA) *PFS*SIN(PSI)**2+(-2*COS (BETA) *SIN(BE 

; TA) *COS (DELTA) **2*PFS*C0S (PSI) +( 1-2*C0S(BETA)**2)*C0S(DELTA)**2 

< * 0MG-2*C0S( BETA ) *COS (DELTA )* *2* 'DIFF (DELTA ,T,l))*SIN(PSI)-SIN(B 

= ETA ) *C0S (DELTA )*SIN (DELTA )*PFS*COS( PSI)* *2+ (-COS (BETA) *COS(DELT 

> A) *SIN (DELTA) *0MG-2*C0S (DELTA) *SIN (DELTA) * 'DIFF(DELTA , T , 1 ) ) *COS 

? (PSI) )*QFS+( COS (DELTA )**2*PFS*C0S (PSI ) -COS (BETA) *COS (DELTA) *SIN 
a (DELTA) *PFS*SIN (PSI) )* 'DIFF (PSI , T , l)+(-SIN (BETA) *C0S (DELTA) *S IN 

1 (DELTA) *PFS**2*C0S (PSI) +C0S (DELTA) **2* 'DIFF(PFS ,T , 1 )+(-COS(BETA 

2 )*COS (DELTA) *SIN (DELTA) *0MG-2*C0S (DELTA) *SIN(DELTA)*'DIFF(DELTA 

3 , T , 1 ) ) *PFS) *SIN (PSI) +COS( BETA) *S IN (BETA) *COS (DELTA) **2*PFS**2*C 

4 OS (PSI ) **2+ (COS (BETA) *COS( DELTA) *SIN (DELTA)* 'DIFF(PFS ,T , 1 )+( (2* 

5 COS (BETA) **2- 1 ) *C0S(DELTA) **2*OMG+2*COS(BETA)*COS (DELTA) **2*'DI 

6 FF (DELTA ,T , 1) )*PFS)*COS(PSI)-SIN(BETA) *COS (DELTA) *SIN(DELTA) * 'D 

7 IFF (OMG , T , 1 )-COS (BETA) *SIN ( BETA )*COS( DELTA) **2*0MG**2-2*SIN (BET 

8 A) *COS (DELTA ) **2* ' DIFF (DELTA ,T , 1 ) *0MG+2* ’ DIFF (BETA ,T , 1 ) *C0S (DEL 

9 TA)*SIN (DELTA)* ' DIFF (DELTA ,T, 1) ) *SM+SIN(BETA) *DELE*SIN (DELTA) *F 

: M* 'DIFF(RFS , T , 1 ) + ( (-SIN (BETA) *DELE*EBETA-COS (BETA) *SIN (BETA) *DE 

; LE**2)*M+( -SIN (BETA) *COS (DELTA )*EBETA-2*C0S( BETA) *SIN(BETA) *DEL 

< E*COS (DELTA) ) *FM)*RFS**2+ ( ( ( (COS (BETA) *DELE*EBETA+ (2*C0S (BETA) * 

= *2-1 ) *DELE**2) *M+ (COS (BETA) *COS (DELTA) *EBETA+ (4*C0S (BETA) **2-2) 

> *DELE* COS (DELTA) ) *FM) *SIN( PSI) +COS (BETA) *DELE*SIN (DELTA) *FM*COS 
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? (PSI))*QFS+C0S(BETA)*DELE*SIN(DELTA)*FM*PFS*SIN(PSI)+(((1-2*C0S 
® ( BETA )**2)*DELE**2-C0S( BETA )*DELE*EBETA)*M+( (2-4*C0S(BETA)**2) * 

1 DELE+COS (DELTA ) -COS ( BETA) *COS (DELTA )*EBETA)*FM)*PFS*COS( PSI ) + ( ( 

2 2*SIN (BETA) *DELE*EBETA+2*COS (BETA) *SIN (BETA) *DELE**2) *M+ (2*SIN ( 

3 BETA) *COS (DELTA ) *EBETA+4*C0S (BETA) *SIN (BETA) *DELE*COS (DELTA) ) *F 

4 M)*0MG+2*SIN ( BETA) *DELE*COS (DELTA)* ' DIFF (DELTA , T , 1 ) *FM)*RFS+ ( ( ( 

5 COS ( BETA ) *DELE*EBETA+DELE**2 )*M+(COS (BETA) *COS (DELTA ) *EBETA+2*D 

6 ELE*COS (DELTA ) ) *FM) *COS (PSI ) -COS (BETA) *DELE*SIN (DELTA) *FM*SIN (P 

7 SI) )* ; DIFF(QFS ,T , l)+( ( COS ( BETA )*S IN ( BETA) *DELE**2*M+2*C0S (BETA) 

8 *SIN (BETA) *DELE*COS (DELTA )*FM) *SIN(PSI)**2+SIN (BETA)*DELE*SIN(D 

9 ELTA) *FM*COS (PSI) *SIN (PSI)+ (-SIN (BETA) *DELE*EBETA*M-SIN (BETA) *C 

: 0S(DELTA)*EBETA*FM)*C0S(PSI)**2)*QFS**2+( ( (( -COS ( BETA) *DELE*EBE 

; TA-DELE* *2 ) *M+( -COS (BETA) *COS (DELTA) *EBETA-2*DELE*C0S (DELTA) ) *F 

< M)*SIN (PSI) -COS (BETA) *DELE*SIN (DELTA) *FM* COS (PSI) ) * ’DIFF (PSI ,T , 
= l)+SIN(BETA)*DELE*SIN(DELTA)*FM*PFS*SIN(PSI)**2+(((-2*SIN(BETA) 

> *DELE*EBETA-2*C0S ( BETA ) *SIN ( BETA) *DELE**2 ) *M+ ( -2*SIN (BETA ) *COS ( 

? DELTA) *EBETA-4*COS (BETA) *SIN (BETA) *DELE*COS(DELTA) ) *FM)*PFS*COS 

® (PSI ) + ( ( ( 1-2*C0S (BETA) **2) *DELE**2-COS (BETA) *DELE*EBETA) *M+ ( (2- 

1 4*C0S ( BETA) **2) *DELE*COS (DELTA) -COS (BETA ) *COS (DELTA ) *EBETA)*FM) 

2 *OMG-2*COS( BETA )*DELE*COS (DELTA)* ’DIFF (DELTA , T, 1)*FM)*SIN (PSI)~ 

3 SIN ( BETA ) *DELE*SIN (DELTA) *FM*PFS*COS (PSI ) **2+ (-COS ( BETA ) *DELE*S 

4 IN (DELTA) *FM*OMG-2*DELE*SIN (DELTA)* 'DIFF (DELTA ,T, l)*FM)*COS(PSI 

5 ) ) *QFS+ ( ( (COS (BETA) *DELE*EBETA+DELE**2) *M+ (COS (BETA) *COS (DELTA) 

6 *EBETA+2*DELE*C0S (DELTA) ) *FM) *PFS*COS (PSI ) -COS (BETA ) *DELE*SIN (D 

7 ELTA) *FM*PFS*SIN (PSI) ) * ’DIFF (PSI , T , 1 ) + (-SIN (BETA) *DELE*EBETA*M- 

8 SIN(BETA)*C0S(DELTA)*EBETA*FM)*PFS**2*SIN(PSI)**2+(-SIN(BETA)*D 

9 ELE*SIN (DELTA) *FM*PFS**2*C0S (PSI ) + ( (COS ( BETA) *DELE*EBETA+DELE** 

: 2)*M+( COS ( BETA ) *COS (DELTA ) *EBETA+2*DELE*C0S (DELTA ) ) *FM) * ' DIFF (P 

; FS,T, l) + ( -COS (BETA) *DELE*SIN (DELTA) *FM*0MG-2*DELE*SIN (DELTA) * ’ D 

< IFF (DELTA , T , 1 ) *FM) *PFS+AHFSY*SIN ( BETA) *DELE*M+AHFSY*SIN (BETA) *C 

= OS (DELTA) *FM) *SIN (PSI )+ (COS (BETA) *SIN( BETA) *DELE**2*M+2*C0S (BET 

> A) *SIN (BETA) *DELE*COS (DELTA) *FM) *PFS**2*C0S(PSI ) **2+ (COS (BETA) * 

? DELE*SIN (DELTA)*FM* J DIFF(PFS ,T, l) + ( ( (COS (BETA) *DELE*EBETA+(2*C0 
fi S (BETA ) **2-1 ) *DELE**2)*M+( COS ( BETA) *C0S (DELTA )*EBETA+(4*C0S( BET 

1 A) **2-2) *DELE*COS (DELTA) ) *FM) *0MG+2*C0S (BETA) *DELE*COS (DELTA) * » 

2 DIFF(DELTA ,T , l) *FM)*PFS-AHFSX*SIN (BETA)*DELE*M-AHFSX*SIN (BETA)* 

3 COS (DELTA) *FM) *COS (PSI ) -SIN ( BETA) *DELE*SIN (DELTA) *FM* J DIFF ( OMG , 

4 T , 1 ) + ( (-SIN (BETA ) *DELE*EBETA-COS (BETA) *SIN (BETA) *DELE**2) *M+ (-S 

5 IN ( BETA) *COS (DELTA) *EBETA-2*C0S (BETA) *SIN (BETA)*DELE*COS (DELTA) 

6 )*FM) *0MG**2-2*SIN(BETA)*DELE*C0S (DELTA) * J DIFF (DELTA ,T , l) *FM*OM 

7 G-MRF2X-MAF2X+AHFSZ*C0S(BETA)*DELE*M+(2*’DIFF(BETA,T,l)*DELE*SI 

8 N (DELTA) * » DIFF (DELTA , T , 1) +AHFSZ*COS (BETA) *C0S (DELTA) ) *FM) /(COS ( 

9 DELTA) **2*SM+DELE**2*M+2*DELE*C0S (DELTA) *FM) 


35 




Appendix B 
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:f <31 INTENTIONALLY BLANK 



SUBROUTINE ORG : 


This subroutine is an interactive routine for simplifying a given expression 
"A" . The call is of the form: 

ORG (A, [STATUS] ) 

where "A" is the expression to be siirplified. After the routine is envoked 
it changes the regular MACSYMA proirpt to ORG> signifying that the program 
is awaiting the expression to be factored out . After the user enters this 
expression, followed by a as usual, the subroutine factors out the 

given expression, gives the part that was factored out (calls it THIS TIME:) . 
Entering a new expression at this point will continue the process by factoring 
out this newly entered expression from the REMAINS part. Entering a "QUIT" 
will quit the routine. The result of the factoring operation is the sum of the 
two parts RUNSUM and RUNREM. If complete factorization is accomplished, RUNREM 
will be zero. 

***** 

*/ 

ORG (EXP, STATUS) : =BLOCK ( [FACT, RUNPART] , 

RUNREM: EXP, 

IF STATUS -NEW THEN RUNSUM: 0, 

DO 

(FACT : READ ( n 0RG> "), 

IF FACT # QUIT THEN 

(if atom (FACT) then 

( if ratcoef (runrem, FACT, 1)=0 then 
runpart : 0 
else 

( runpart : FACT*TRIGS IMP < 

RATSIMP (ratcoef (runrem, FACT, 1) ) ) , 
runrem: ratcoef (runrem, FACT, 0) , 
runsum: runsum+ runpart) ) 

else 

(RUNREM :RATSUBST( SUB, FACT, RUNREM) , 

IF RATCOEF (RUNREM, SUB, 1)=0 THEN 
(RUNPART: 0, 

runrem: subst (FACT, sub, runrem) ) 

ELSE 

(RUNPART : SUB*TRIGSIMP ( 

RATSIMP (RATCOEF (RUNREM, SUB, 1) ) ) , 

RUNREM : RATEXP AND (RATCOEF (RUNREM, SUB, 0) ) , 

RUNREM: SUBST (FACT, SUB, RUNREM) , 

RUNPART : SUBST (FACT, SUB, RUNPART) , 

RUNSUM : RUNSUM+RUNPART ) ) , 

PRINT(”THIS TIME: ", RUNPART)) 
else 

RETURN ("END ORGANIZE"))); 
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Table 1. Aircraft inertias (slug x ft 2 ) 


I XX 

6852 

lyy 

37961 

Izz 

38260 

Ixz 

4379 


Table 2. Aircraft configuration data for hover maneuvers 


Maneuver 

Main input 

Airspeed 

(kts) 

Weight 

(ib) 

C.G. 

(in.) 

Altitude 

(&) 

Start time 
(sec) 

4020 

Lateral 

Hover 

15399 


mmrm® 

2.8 

4017 

Long. 

Hover 

15475 



1.0 

4014 

Dir. 

Hover 

15619 



2.0 

4028 

Heave 

Hover 

15184 



2.8 


Table 3. Aircraft configuration data for 60 kts maneuvers 


Maneuver 

Main input 

Airspeed 

(kts) 

Weight 

(ib) 

C.G. 

(in.) 

Altitude 

(ft) 

Start time 
(sec) 

4927 

Lateral 

71 

14957 

204.80 

6000 

0.0 

4912 

Long. 

72 

15580 

204.69 

6000 

0.0 

5105 

Dir. 

70 

15811 

204.66 

6000 

0.0 

5112 

Heave 

75 

15372 

204.74 

6000 

0.0 


Table 4. Aircraft configuration data for 120 kts maneuvers 


Maneuver 

Main input 

Airspeed 

(kts) 

Weight 

(ib) 

C.G. 

(in.) 

Altitude 

(ft) 

Start time 
(sec) 

4317 

Lateral 

128 

15246 

204.65 

6000 

0.0 

4310 

Long. 

130 

15887 

204.55 

6000 

0.0 

4320 

Dir. 

135 

15103 

204.67 

6000 

0.0 

4326 

Heave 

131 

14740 

204.73 

6000 

0.0 
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Table 5. (a) Stability derivatives in hover 



Flight 

BEMAP 

FLYRT 

Derivative 

Par am. value 

Param. value 

Param. value 


-0.02000 f 

-0.03290 

-0.03320 

x w 

0.000 + 

0.02280 

0.02790 


0.000 + 

1.715 

3.023 

x v 

0.000 + 

-7.000E-03 

0.01210 

x p 

0.6849 

-1.117 

-0.4413 

x r 

0.000 + 

-0.02830 

-0.03220 

Z u 

0.000 + 

6.300E-03 

4.500E-03 

Zyj 

-0.1220 

-0.3171 

-0.2921 

Z q 

0.000 + 

0.3261 

0.1199 

z v 

0.000 + 

1.100E-03 

-3.900E-03 

z v 

0.000 + 

-0.1442 

-0.2241 

z r 

0.000 + 

2.357 

2.171 

M u 

8.443E-04 

-4.500E-03 

-3.700E-03 

Myj 

-5.142E-03 

1.000E-03 

-6.000E-04 

Mg 

-0.4192 

-0.6069 

-0.5521 

M v 

7.103E-03 

2.800E-03 

-1.400E-03 

M p 

-0.2272 

0.2142 

0.08750 

Mr 

-0.09000 f 

0.01510 

0.01160 

Yu 

0.000 + 

0.01790 

-8.300E-03 

Yw 

0.000 + 

-1.900E-03 

-0.01740 

Y, 

0.000 + 

-1.375 

-1.203 

Yy 

-0.2788 

-0.05790 

-0.05020 

Yp 

-1.560 

-2.672 

-2.938 

Y r 

0.000 + 

0.6762 

0.4704 

Lu 

0.000 + 

0.02360 

-6.900E-03 

Lw 

0.000 + 

-2.100E-03 

-0.01150 

L q 

1.040 f 

-1.171 

-0.9598 

Ly 

-4.247E-03 

-0.01820 

-0.01710 

L p 

-1.828 

-2.977 

-2.424 

l'jj' 

0.000 + 

0.4681 

-0.07800 

N u 

0.000 + 

2.700E-03 

-1.800E-03 

N w 

0.000 + 

2.000E-04 

-7.000E-04 

w. 

0.000 + 

-0.1570 

-0.1114 

Ny 

3.008E-03 

6.400E-03 

6.600E-03 

Np 

-0.3085 

-0.2605 

-0.3376 

N r 

-0.2702 

-0.2291 

-0.3176 


-T Eliminated during model structure determination. 
•[Fixed value in model. 
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Table 5. (b) Control derivatives in hover 



Flight 

BEMAP 

FLYRT 

Derivative 

Par am. value 

Param. value 

Param. value 

Xlon 

-1.483 . 

-1.519 

-1.841 

Xiat 

-0.1939 

0.04810 

-0.07600 

Xped 

0.000 + 

-0.07430 

-0.06320 

Xcol 

0.8355 

0.6636 

0.6897 

Zion 

0.000 + 

-0.09910 

-0.1563 

Zlat 

0.000 + 

3.400E-03 

-0.03070 

Zped 

0.000 + 

-8.000E-03 

-1.000E-04 

Zcol 

-5.430 

-7.306 

-7.863 

Ml on 

0.2353 

0.2625 

0.3144 

M lat 

0.05917 

-0.02710 

0.01330 

Mp e d 

0.000 + 

0.01750 

0.01810 

M co l 

0.000 + 

5.000E-04 

1.500E-03 

Ylon 

0.000 + 

-0.2022 

-0.1123 

Ylat 

0.4958 

1.266 

1.287 

Xped 

-2.791 

-1.796 

-1.708 

Ycol 

-0.8557 

-0.2165 

-0.1659 

^lon 

-0.1042 f 

-0.1370 

-0.1000 

Llat 

0.8341 

1.184 

1.145 

Lped 

-0.4008 

-0.5612 

-0.3560 

Lcol 

0.000 + 

-0.07360 

-0.02640 

Nlon 

0.000 + 

-0.01770 

-0.01310 

Niat 

0.1041 

0.1553 

0.1501 

Nped 

0.4943 

0.5763 

0.5887 

Xcol 

0.2662 

0.1941 

0.2915 

T lon 

0.08830 

0.08830 

0.08830 

Tlat 

0.1206 

0.1206 

0.1206 

Tped 

0.07951 

0.07951 

0.07951 

Tool 

0.06107 

0.06107 

0.06107 


+ Eliminated during model structure determination. 
fFixed value in model. 
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Table 6. (a) Stability derivatives at 60 kts 



BEMAP 

FLYRT 

Derivative 

Par am. value 

Par am. value 


-0.03390 

-0.03530 


0.02470 

0.01900 

x q 

-0.03360 

1.023 

x v 

0.02040 

0.05380 

x p 

-1.013 

-0.1456 

x r 

-0.3169 

-0.1559 

Z u 

-0.06070 

-0.05480 

Z w 

-0.5419 

-0.5034 

z q 

98.40 

101.3 

Z v 

-0.4638 

-0.4439 

Zp 

4.714 

-0.5447 

Z r 

2.557 

1.797 

M u 

1.700E-03 

1.500E-03 

M w 

2.300E-03 

1.500E-03 

M q 

-1.368 

-0.8738 

M v 

3.300E-03 

-2.200E-03 

M p 

0.1127 

0.04590 

Mr 

0.02980 

0.01120 

Yu 

5.400E-03 

5.700E-03 

Y w 

1.200E-03 

-0.03760 

Y, 

-1.087 

* -1.031 

Y v 

-0.08810 

-0.1706 

Yp 

1.615 

-0.6978 

Y r 

-100.3 

-100.3 

Lu 

-2.700E-03 

3.700E-03 

L w 

0.01720 

-0.01790 

L q 

-0.8584 

-1.181 

L v 

-9.700E-03 

-0.08480 

Lp 

-4.290 

-3.067 

J.jj' 

-0.1564 

-0.1169 

N u 

-3.600E-03 

-3.600E-03 

N w 

-2.000E-03 

-7.800E-03 

N q 

0.1280 

-0.02520 

N v 

6.300E-03 

-4.800E-03 

N p 

-0.5311 

-0.4812 

N r 

-0.4518 

-0.5120 
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Table 6. (b) Control derivatives at 60 kts 



BEMAP 

FLYRT 

Derivative 

Param. value 

Param. value 

Xlon 

-1.615 

-1.733 

Xiat 

-0.1791 

-0.06850 

Xped 

0.1434 

0.1150 

Xcol 

0.1470 

0.1416 

Zion 

-1.630 

-2.462 

Zlat 

0.1510 

-0.07630 

Zped 

1.200E-03 

8.000E-04 

Zcol 

-6.982 

-6.690 

Mlon 

0.3700 

0.3316 

Mlat 

0.03360 

0.01500 

Mped 

-0.01550 

-5.400E-03 

M col 

0.09080 

0.07190 

Ylon 

-0.1311 

-0.3392 

Ylat 

0.9025 

1.227 

Yped 

-1.163 

-1.254 

Ycol 

0.04740 

-0.1801 

Lion 

-0.08710 

-0.2402 

Llat 

1.263 

1.100 

Lped 

-0.2921 

-0.2558 

LqoI 

0.1696 

-0.1233 

Xlon 

-0.02670 

-0.08140 

Xiat 

0.1495 

0.1482 

Xp e d 

0.3740 

0.4367 

Xcol 

0.09920 

0.1018 

T lon 

0.08830 

0.08830 

Tlat 

0.1206 

0.1206 

Tped 

0.07951 

0.07951 

Tool 

0.06107 

0.06107 
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Table 7. (a) Stability derivatives at 120 kts 



BEMAP 

FLYRT 

Derivative 

Param. value 

Param. value 


-0.05680 

-0.05520 

A™ 

0.02640 

0.01360 

X q 

10.37 

10.29 

x v 

0.04610 

0.05670 

x p 

-0.9485 

-0.5177 

x r 

-0.06490 

-0.05550 

Z u 

0.01620 

0.02210 

Z w 

-0.7376 

-0.6446 

Z, 

201.0 

201.8 

z v 

-1.239 

-1.138 

Zp 

1.775 

-2.334 

z r 

3.022 

2.229 

M u 

2.200E-03 

1.800E-03 

Myj 

2.600E-03 

2.500E-03 

M q 

-1.430 

-1.077 

M v 

7.400E-03 

4.400E-03 

M p 

0.1902 

0.1295 

M r 

0.02380 

0.02190 

Y 

1 u 

0.01420 

0.01390 

y 

1 w 

-3.800E-03 

-0.03580 

y. 

-1.044 

-0.6811 

Y v 

-0.1422 

-0.2552 

Yv 

-9.555 

-10.40 

Y r 

-200.7 

-200.8 

L u 

-4.000E-04 

3.200E-03 

Lyj 

0.01440 

-9.500E-03 

L„ 

-1.129 

-0.8714 

L v 

-3.600E-03 

-0.09040 

Lp 

-3.626 

-3.168 

Lj* 

-0.07620 

-0.06570 

N u 

-5.200E-03 

-5.900E-03 

N w 

-3.000E-04 

-2.500E-03 

N g 

0.06130 

2.300E-03 

N v 

8.300E-03 

2.000E-03 

N p 

-0.2875 

-0.2754 

N r 

-0.6668 

-0.7556 
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Table 7. (b) Control derivatives at 120 kts 



BEMAP 

FLYRT 

Derivative 

Par am. value 

Param. value 

Xion 

-1.479 

-1.771 

Xiat 

-0.1824 

-0.06870 

Xped 

0.09060 

0.08340 

X co l 

7.200E-03 

-0.04460 

Zion 

-5.472 

-5.541 

Zlat 

0.07900 

-0.06750 

Zped 

0.03120 

0.01020 

Zcol 

-8.921 

-7.812 

Mlon 

0.4585 

0.3965 

M[ at 

0.02830 

0.01510 

Mped 

-0.01170 

-1.800E-03 

Mcol 

0.1882 

0.1579 

Ylon 

-0.08940 

-0.4851 

Yl at 

0.9004 

1.267 

Yped 

-1.484 

-1.569 

Y coi 

-7.400E-03 

-0.3069 

Zion 

0.05590 

-0.2438 

Zlat 

1.280 

1.111 

Zped 

-0.3844 

-0.3294 

Zcol 

0.2402 

-0.05930 

Nlon 

-0.01900 

-0.07390 

Nlat 

0.1541 

0.1491 

Nped 

0.4676 

0.5437 

Xcol 

0.1402 

0.1629 

r lon 

0.08830 

0.08830 

T lat 

0.1206 

0.1206 

7~ped 

0.07951 

0.07951 

Tool 

0.06107 

0.06107 
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Figure 1 . AH-64A three view illustration. 



(a) *2 

Figure 2. (a) Coordinate systems used in rotor derivations. 


49 



2 b 



Vrs 


systems (drawing not to scale), (c) transformation from fixed to 

rotating-shaft coordinates. 
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o Flight data 





(d) Airspeed (knots) 


o Flight data 
BEMAP 



o Flight data 




(g) Airspeed (knots) 


Figure 3. (a) Variation of trim pitch attitude with airspeed, (b) variation of trim collective with airspeed, 
(c) variation of trim longitudinal cyclic with airspeed, (d) variation of trim roll attitude with airspeed, 

(e) variation of trim lateral cyclic with airspeed, (f) variation of trim pedals with airspeed, (g) variation 
of trim horizontal stabilator incidence with airspeed. 
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Rotor speed (rpm) Yaw rate (deg/sec) Pitch rate (deg/sec) Roll rate (deg/sec) 


30.0 


BEMAP V-1.0 

FLYRT V-5.13 





(b) Time (sec) 

Figure 4. (b) Continued. 
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Horizontal stabiiator Normal acceleration (g) Lateral velocity (ft/sec) Longitudinal velocity (ft/sec) 

incidence (deg) 






(c) Time (sec) 


Figure 4. (c) Concluded. 
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Collective (%) Pedal (%) Longitudinal cyclic (%) Lateral cyclic (%) 
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BEMAP V-1.0 
FLYRT V-5.13 
Trends flight - test data 
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Horizontal stabilator 

incidence (deg) Normal acceleration (g) Lateral velocity (ft/sec) Longitudinal velocity (ft/sec) 
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Figure 6. (a) Responses to directional doublet in hover. 
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Rotor speed (rpm) Yaw rate (deg/sec) Pitch rate (deg/sec) Roll rate (deg/sec) 






Figure 6. (b) Continued. 
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Horizontal stabilator 

incidence (deg) Normal acceleration (g) Lateral velocity (ft/sec) Longitudinal velocity (ft/sec) 


BEMAP V-1.0 
FLYRT V-5.13 
Trends flight - test data 



Figure 6. (c) Concluded. 
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BEMAP V-1.0 
FLYRT V-5.13 
Trends flight - test data 





















Time (sec) 


Figure 9. (a) Responses to longitudinal doublet at 60 kts 
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Figure 10. (a) Responses to directional doublet at 60 kts. 
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Figure 1 1 . (b) Continued. 











Figure 12. (a) Responses to lateral doublet at 120 kts. 
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Yaw rate 



Figure 12. (b) Continued. 
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Collective {%) Pedal {%) Longitudinal cyclic (%) Lateral cyclic (%) 
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BEMAP V-1.0 
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Figure 13. (a) Responses to longitudinal doublet at 120 kts. 
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Rotor speed (rpm) Yaw rate (deg/sec) Pitch rate (deg/sec) Roll rate (deg/sec) 






(b) Time (sec) 

Figure 13. (b) Continued. 
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Horizontal stabilator 

incidence (deg) Normal acceleration (g) Lateral velocity (ft/sec) Longitudinal velocity (ft/sec) 
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10 


BEMAPV-1.0 
FLYRT V-5.13 
Trends flight - test data 
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Figure 15. (b) Continued 
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BEMAP, input at swashplate 

BEMAP, input at stick 

T rends flight - test data 





Time (sec) 

Figure 16. Lateral doublet at 60 kts, input at stick. 
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Figure 17. (a) Comparison of BEMAP responses with flight at hover. 
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Figure 17. (b) Continued. 
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Flight results 
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Figure 18. (a) Comparison of BEMAP responses with flight at 60 kts. 
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Flight results 

BEMAP V-1.0 




Figure 18. (b) Continued. 
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Flight results 

BEMAP V-1.0 




Figure 1 8. (c) Concluded. 


96 






97 


REPORT DOCUMENTATION PAGE 


Form Approved 
OMBNo. 0704-0188 


Public reporting burden for this collection of information is estimated to average 1 hour per response, including the time for reviewing instructions, searching existing data sources, 
gathering and maintaining the data needed, and completing and reviewing the collection of information. Send comments regarding this burden estimate or any other aspect of this 
collection of information, including suggestions for reducing this burden, to Washington Headquarters Services, Directorate for information Operations and Reports, 1215 Jefferson 
Davis Highway, Suite 1204, Arlington, VA 22202-4302, and to the Office of Management and Budget, Paperwork Reduction Project (0704-0188), Washington, DC 20503. 


1. AGENCY USE ONLY (Leave blank) I 2. REPORT DATE 


April 1995 


3. REPORT TYPE AND DATES COVERED 

Technical Memorandum 


4. TITLE AND SUBTITLE 


Development and Validation of a Blade-Element Mathematical 
Model for the AH- 64 A Apache Helicopter 


6. AUTHOR(S) 


M. Hossein Mansur 


5. FUNDING NUMBERS 


505-59-36 


7. PERFORMING ORGANIZATION NAME(S) AND ADDRESS(ES) 

Aeroflightdynamics Directorate, U.S. Army Aviation and Troop 
Command, Ames Research Center, Moffett Field, California 


8. PERFORMING ORGANIZATION 
REPORT NUMBER 


A-95034 


9. SPONSORING/MONITORING AGENCY NAME(S) AND ADDRESS(ES) 

National Aeronautics and Space Administration 

Washington, DC 20546-0001 and U.S. Army Aviation and Troop 

Command, St. Louis, MO 63120-1798 


10. SPONSORING/MONITORING 
AGENCY REPORT NUMBER 

NASA TM- 108863 
USAATCOM TR-94-A-022 


11. SUPPLEMENTARY NOTES 



Point of Contact: M. Hossein Mansur, Ames Research Center, MS 211-2, 
Moffett Field, CA 94035-1000; (415) 604-6037 


12a. DISTRIBUTION/AVAILABILITY STATEMENT 

Unclassified-Unlimited 
Subject Category - 01 


13. ABSTRACT (Maximum 200 words) 

A high-fidelity blade-element mathematical model for the AH-64A Apache Advanced Attack 
Helicopter has been developed by the Aeroflightdynamics Directorate of the U.S. Army’s Aviation and 
Troop Command (ATCOM) at Ames Research Center. The model is based on the McDonnell Douglas 
Helicopter Systems’ (MDHS) Fly Real Time (FLYRT) model of the AH-64A (acquired under con- 
tract) which was modified in-house and augmented with a blade-element-type main-rotor module. This 
report describes, in detail, the development of the rotor module, and presents some results of an exten- 
sive validation effort. 


14. SUBJECT TERMS 


Blade-element; Helicopter, Mathematical model 


15. NUMBER OF PAGES 

98 


16. PRICE CODE 

A05 


17. SECURITY CLASSIFICATION 118. SECURITY CLASSIFICATION 119. SECURITY CLASSIFICATION 120. LIMITATION OF ABSTRACT 

OF REPORT I OF THIS PAGE I OF ABSTRACT 1 


Unclassified 


NSN 7540-01-280-5500 


Unclassified 


Standard Form 298 (Rev. 2-89) 

Prescribed by ANSI Std. Z39-18 






















